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SUMMARY 


The thesis deals with the development of an analytical technique 
for predicting the acoustic performance of turbofan inlets carrying a 
subsonic axi syiwnetr ic steady floi. The Finite Elet^ent Method (FEM) in 
combination with the Method of Weighted Residuals has been chosen as 
the solution technique for predicting the acoustic properties of vari- 
able area, annular ducts with or without acoustic treatments along 
the i r walls. 

An approximate solution for the steady inviscid flow field is 
obtained using an integral method for calculating the incompressible 
potential flow field in the inlet with a correction to account for com- 
pressibility effects. 

The accuracy of the finite elwent technique in predicting the 
acoustic properties of annular ducts has been checked by comparison with 
available analytical solutions for the problems of plane and spinning 
wave propagation through a hard walled annular cylinder with a cons;.) it 
mean flow. 

For a fixed number of triangular elements in the finite element 
scheme, the number of nodes per wavelength decreases as the frequency 
of oscillation Increases, resulting In a loss in accuracy in the numeri- 
cal results. Accuracy at higher frequency can be recovered by increas- 
ing the number of nodes per wavelength, either by increasing the number 
of linear elements (3"nodes per triangle) or by utilizing a ftx^re 
elaborate description of the dependent variable over each element as is 



XV 


the case when quadratic el^nents (6-nodes per triangle) are used. It 
Is clearly demonstrated for the case of plane wave propagation in a 
hardwalled annular cylinder containing a uniform steady flow that for 
an equal number of triangular elements, quadratic representation is 
superior to linear representation In handling high frequency wave propa- 
gation the reason being, the use of quadratic elenents effectively 
doubles the number of nodes per wavelength as compared to an equal num- 
ber of linear elements. 

The accuracy of the developed finite element solution is again 
confirmed by compar ing . i ts predictions with those obtained by a finite 
difference solution approach. In this study, the results of the FEM 
for the case of sound propagation in an acoustically lined annular 
cylinder of a large radius of curvature and thin annular spacing carry- 
ing a uniform steady flow are compared with the comprehensive results 
of Baumeister for the equivalent case of a rectangular duct carrying 
a uniform flow. The duct attenuation and the acoustic pressure dis- 
tributions obtained by both the methods are found to be generally in 
good agreement. 

Since exact values of the impedance at the open end of a duct 
are generally not known, simple impedance boundary conditions such as 
"no reflection" impedance conditions are used to date in most cases. 
However, in the absence of a steady flow an integral solution approach 
developed by Bell, Meyer and Zinn provides the "exact" impedance con- 
dition at the open end of the duct. In order to ccxnpare the solution 
technique for the study of non-un I form ducts without steady flows, the 
finite element solutions were required to satisfy the Impedance boundary 



XV i 

conditions obtained using the above mentioned integral solution approach 
and the resulting finite element solutions for the QCSEE (Quiet, Clear, 
Short-haul, Experimental, Engine) inlet are coinpa red wi th those obtained 
for the same inlet, using the integral solution approach. The results 
obtained by the finite elanent program using the "exact" Impedance con- 
dition are found to be in excellent agreement with the results of the 
integral solution approach for plane wave propagation in hard and soft 
walled QCSEE inlet carrying no mean flow. The wave structure within 
the inlet obtained by prescribing the simple "no reflection" impedance 
boundary condition is quite different from the one obtained by prescrib- 
*''1 the "exact" impedance boundary condition for low frequencies of wave 
P -'Psgation. However, as frequency increases the two Impedance boundary 
conditions approach one another and so do the corresponding wave struc- 
tures. The propagation of high frequency plane waves in the inlet could 
be approximated as the wave propagation in a ray tube where the reflected 
component is negligible, in which case, the "no reflection" impedance 
condition indeed approaches the "exact" impedance condition and so also 
the corresponding wave structures. 

Results are presented comparing low frequency plane wave propa- 
gation through the hard walled QCSEE inlet containing a one-dimensional 
steady flow with the same inlet containing a fully two-dimensional 
ax i symmetr i c steady flow. It is shown that when one-dimensional steady 
flow is assumed to exist In the inlet, the plane wave propagates with 
relatively little distortion. However, propagation of a plane wave 
through the fully two-dimensional flow field in the inlet produces 
severe distortions due to the excitation of higher order modes. 
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To estimate the sensitivity of Inlet curvature, center body and 
mean flow gradients on duct attenuation, acoustic calculations for soft 
walled QCSEE Inlet, a straight cylinder and a Bellmouth Inlet of the 
same overall dimensions are performed for various frequencies for the 
cases of zero mean flow and fully two-dimensional axisymmetric mean flow. 
The liner Impedance values chosen are the same for each inlet at a 
particular frequency and they correspond to the near optl.um impedance 
values for a plane pressure wave Input into an infinitely long circular 
duct with zero mean flow. At high frequencies focusing of the sound 
wave toward the duet axis occurs for both zero mean flow and non-zero 
mean flow cases and hence the duct attenuation falls rapidly with 
increase in frequency for all the three inlets. Attenuation of lov 
frequency plane waves is found to be more sensitive to inlet curvature, 
center body and mean flow gradients (for the non-zero mean flow case) 
compared to that of high frequency plane waves. The approximate method 
proposed by Rice to obtain the near optimum liner impedance values for 
the flow case is found to be reasonable for high frequency plane waves 
but leads to a rapid drift in the optimum values for the low frequency 
plane waves resulting in a considerable reduction of duct attenuation. 
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CHAPTER I 
INTRODUCTION 
I . 1 General 

As evidenced by the recent surge in conmercial aviation mainly due 
to Federal deregulation of price control on air travel by the US Civil 
Aeronautics Board (C.A.B.)* the aircraft industry has come up with a 
new generation of aircrafts to meet the record high demands of domestic 
and international air travel. The new series Boeing 767 and 757 air- 
crafts, the inodified Lockheed L 1011 aircraft and the shorter version 
of French Airbus A300B aircraft will be put Into airline service in 
early I980s to satisfy the growing needs of air transportation. With 
the increase in air traffic and consequent congestion in airport opera- 
tions, it has become increasingly apparent to government, industry and 
other research and development organizations that major efforts must be 
undertaken to improve the general community environment affected by 
the convnercial and also military air transportation. The recent legis- 
lation of the Federal Act FAR 36 (1978) calls for very stringent meas- 
ures to be taken by aircraft industries to reduce the overall noise 
levels of their aircraft. Since the aircraft engines are the major 
contributors to the overall aircraft noise, one has to develop a capa- 
bility to accurately predict the noise levels due to various components 
of the aircraft engines, which in turn will enable the engine designer 
to develop a viable design capable of meeting the new lower noise levels 
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without adversely affecting the performance of the power plant. 

Aircraft generated noise sources can be divided Into two groups; 
namely the externally generated noise produced by the Jet exhaust, and 
the Internally generated noise that Is primarily due to the rotating 
turbomachinery blades and the combustion process. The utilization of 
energy efficient, high bypass ratio turbofan engines In the present 
day civil and military aircraft reduce the jet velocity and hence Jet 
noise compared to earlier low bypass turbofan or turbojet engines. 

The pure tone turbofan noise emitted by these engines now becomes the 
major source of noise pollution. Sound absorbing liners have been 
designed and utilized in the engine inlets to reduce the knitted turbo- 
fan noise mainly on the basis of costly trial and error development 
programs. Hence, the need for an efficient analytical technique which 
can predict the acoustic properties of complex shaped Inlet ducts with 
sound absorbing liners and carrying a multidimensional mean flow for a 
variety of practical noise source inputs is more apparent now than ever 

before. The development of such an analytical technique utilizing the 
f 

Finite Element Method (to be denoted henceforth by FEM) is the object 
of this work. 

1.2 Literature Review 

An extensive survey of the acoustics of aircraft inlets has been 

presented by Nayfeh, Kaiser, and Tellonis^ including a comprehensive 

bibliography. Of particular interest are the use of numerical methods 

2 3 

such as finite differences and integral techniques in the study of 
acoustic propagation in variable area hard walled ducts without mean 
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flow. Due to the additional complications created by the presence of 
a steady flow, most studies of sound propagation in annular duces with 
varying cross sectional area containing a mean flow enploy one or more 

L 

simplifying assumptions such as one- dimensional mean flow, quasi-one 
din«nsional mean flow,^ or slowly varying cross sectional area.^ In 
many instances, however, practical considerations call for the use of 
relatively short ducts having large transverse and streamwise velocity 
gradients. Under such conditions, the predictions of existing theoreti- 
cal approaches in which the mean flow is assumed to be one-dimensional 
or nearly one-dimensional is open to question. Thus, there is a need 
for an analytical method that can determine the acoustics of duct sys- 
tems involving multidimensional flows. 

The relative merits of the various numerical techniques are 

evaluated in this section. The application of the method of finite 

2 

differences by Baumeister, to solve the wave propagation problem is 
restricted to the simple geometry of a rectangular duct carrying a uni- 
form steady flow. The finite difference method employs a uniform rec- 
tangular grid structure which can not adequately represent complex duct 
geometries encountered in aircraft inlet design. Beil, Meyer and Zinn^ 
have developed an integral approach using Green's functions to solve 
the Helmholtz equation for an arbitrarily shaped body. In this method, 
the dimensionality of the problem is reduced by one (e.g. a three dimen- 
sional problem is reduced to a two dimensional surface integral equation, 
an axi-symmetr ic problem is reduced to a line integral equation, etc. 

The integral solution technique is, however, limited to dealing with 
acoustic problems involving either no or constant mean flow. The Galerkin, 
which is a special application of the method of weighted residuals, was 
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applied by Unruh and Eversman^ to solve the wave equation in a duct, 
in this method the solution is sought in terms of a complete set of 
orthogonal basis functions which individually satisfy the natural bound- 
ary conditions and collectively satisfy the forced boundary conditions 
of the problem. The unknown coefficients of the basis functions are 
solved for by requiring the resulting residues to be orthogonal to each 
of the basis functions. This yields a system of linear algebraic equa- 
tions for the coefficients which can be solved by standard matrix methods. 
The disadvantage of the method is that a new set of basis functions has 
to be determined by solving the corresponding eigenvalue problem for 
each frequency. For ducts with complex geometries and acoustic liner 
configurations the evaluation of the basis functions is quite compli- 
cated and time consuming. Hence, a parametric study of ducts for a 
range of frequencies and liner configurations can not be efficiently 
conducted by employing this method. The method of asymptotic expansions^ 
has been used to study the problem of wave propagation in ducts which 
slightly deviate from the simple geometry of a cylinder or a rectangle 
and the resulting mean flow also slightly deviates from a one dimen- 
sional flow. This method obviously can not be used if the duct has 
large variations in shape resulting in a multi dimensional mean flow. 

The finite element method (to be denoted by FEM) which has its 
roots in the field of structural mechanics has been gaining popularity 
in other fields of mechanics because of its ability to deal with complex 
"real life" geometrical shapes and to handle mixed boundary conditions. 

g 

The application of the FEM to the solution of Helmholtz's equation in 
a waveguide indicates the versatility of this method in treating various 
two and three dimensional problems subject to mixed boundary conditions. 

In the finite element method the region of interest is subdivided into 
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a HLNnber of elements of simple shapes (e.g. triangles for t\rtO dimen* 
si<Mial problems and tetrahedrons for three dimensional problems) of ] 

variable sizes which can be arranged to represent complex shapes. This | 

flexibility to utilize a variable mesh distribution of the elements 
for an efficient yet adequate representation of the geometry makes the 
FEM a powerful numerical tool. The FEM is generally thought of in 
connection with variational problems. However, since a variational 
form of the problem of duct acoustics has not been established, the 
governing differential equations are converted to integral equations by 
the Galerkin process. 

1.3 Objectives of Research 

The research work described in this work was initiated for the 
purpose of predicting the acoustic properties of practical turbofan 
inlets carrying high subsonic Mach number mean flows. The inlet config 
urations chosen for study are the Quiet, Clean, Short-haul Experimental 

Q 

Engine (to be henceforth denoted as QCSEE) inlet^ and the Bellmouth 
iniet^*^ (see Figures la and lb). The QCSEE inlet was designed jointly 
by General Electric Co., Douglas Aircraft Company and NASA Lewis 
Research Center for the purpose of developing the technology for a 
turbofan engine intended for application to short take-off and landing 
(STOl) aircraft. The Bellmouth inlet has been used extensively at 
NASA Lewis Research Center as a standard reference base to compare 
experimentally the acoustic performance of realistic aircraft inlets 
because of its simple geometry. The Bellmouth inlet acoustic studies 
are conducted for static cases, the mean flow being created by suction 


t 




Figure Ib. The Bellmouth Inlet Geometry. 
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at the fan plane of the inlet (sea Figure ib). 

To achieve the objectives of this research work, it is necessary 
to develop solution techniques that are capable of predicting the acous* 
tic properties of variable area, annular ducts with or without acoustic 
treatments along their walls, subjected to a variety of practical sound 
excitation conditions. Such a solution technique should be capable of 
properly accounting for the reflection processes at the Inlet entrance 
plane, the space dependence of the noise source at the fan plane (see 
Figure I), the odd shape geometry of the ducts under consideration, the 
mul tidimensionai i ty of the steady flow and the mixed boundary conditions 
(e.g., partial lining) at the duct wails. The above requirements 
obviously preclude the development of an analytical solution for the 
duct under consideration and one must resort to the develojxnent of an 
efficient nurr^rical solution approach. It will be slK}wn in this work 
that the application of the Finite Element Method (FEM) (see Reference 
II for a general discussion of this method), In conjunction with the 
Method of Weighted Residuals can indeed provide the needed solutions. 
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CHAPTER II 

FORMULATION OF THE PROBLEM 
2.1 PlffTcntUl Equations 

To develop the needed solution technique, the problem of acous* 
tic wave propagation through either one of the duct configurations shown 
in Figure I will be considered. The duct carries a two dimensional, 
axisymmetric mean flow which is asswned to be inviscid, ncxi heat con- 
ducting and irrotational . Body forces are neglected. To derive the 
needed nondinwnslonal conservaticm equations, velocities, lengths and 
time are respectively rHjrmalized with respect to the wnbient speed of 
sound c*, a characteristic duct diameter d*, and The density p 

and pressure p are respectively normalized with the ambient density 
and Pq c^. The velocity potential 0 is normalized with respect to c* d* 

and the frequency with respect to c*/d*. Under these conditions. It 
12 

can be shown that the behavior of the flow in the duct is described 
by the following nonlinear partial differential equation for the flow 
potential 0: 

2 

js j 2 2 

2 + fr - c V 4. ( 1 ) 

dt 

where 

■ I ■ (y"I) J 74»*7i|i] 

and Y is the ratio of specific heats. Rewriting Eq. (I) In a cylindrical 


coordinate system (f.e.» r, e* z) with the z*axis coinciding with the 
turbofan iniet axis yields: 


^rr r ^2 ^zz' ^tt 


2 

2 A A 4 — sSpBi + 24(li.-t-4 4 

’r^rt 2 ^z^zt ^r ^rr 

r 

♦ ^ + fWr« + 2* ♦ ♦ 

J* ^z ’zz 2 ^r’^z’^rz 


j 

r r 


( 2 ) 


where the subscripts indicate partial differentiation with respect to 
the subscripted variables. 

To obtain the needed acoustic solutions, the flow potential is 
rewritten as the sum of a steady-axi-symmetr ic mean flow potential <p{r,z) 
and an acoustic potential 0* (r,e,z,t); that is 


0(r,e,z,t) - ^(r,z) ^'(r,e,z,t) (3) 

Because of the rotational nature of the fan and compresser , ^ ^ 
they tend to generate sound that is characterized by spinning acoustic 
modes. In order to account for spinning modes, the acoustic potential 
is assumd to have the following form 

♦ (r.e.z.t) - ♦(r, (I>) 


where 4>(r,z) is a complex quantity; that is 


Substituting Eq. (3) into Eq. (2), neglecting the nwilinear products of 
the acoustic quantities and separating the resulting equation into its 
real and imaginary parts leads to the derivation of the following two 
I inear, coupled, partial differential equations for ^ and h 


and 


* ‘'*'*^*M * ^♦rVrz 
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-2 
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+ [ui “ m “*!f] 0 “ 2u)6 d “ 2u)0 (p 

2 r r 2 2 
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(c^ - + (c^ - - 2i d ♦ 

\w rr 2 ' 22 2 r r 2 

-2 

+ 4-^- (y-I) -f- 

Kk - 

+ 1 -{ y + i (y-i)4>^^^^’ (y-f) —jr^ 

2 2 ' — - 

+ [ui "in ]d ■4' 2<»J^ 0 ^ 2(i)d djj 

r 

+ w(Y-I)[i^^ ■*■ *7?" ° 


where 
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It should be pointed out that In the develofMnent of Eqs. (6) and (7) 
the following expression for the sound velocity perturbation (c ) ‘ has 
been used: 

(c^)' ■ -(y-I) 1-Iti)^' + 

Since the steady flow is axi symmetric, the linearized 9-momentum equa- 
tion can be integrated to give a relationship between the pressure and, 
velocity potential: 

p' -p(-iu^' + (8) 

2.2 Boundary Conditions 

Before proceeding to obtain the needed solutions, the geometry 
and boundary conditions for the problem under consideration must be 
established. The geanetry of typical turbofan inlets are shown in Fig- 
ure I where due to the geometry of the inlets, only a single meridional 
plane Is shown. The boundary of the inlets may be divided into three 
distinct sections, each described b> a different boundary condition. The 
inlet exit plane represents the interface between the inlet and the 
remainder of the engine; it is referred to as the inlet exit plane as it 
is the location where the steady flow leaves the inlet. This plane also 
represents the location at which the fan-compressor noise is Introduced 
into the inlet. In view of the earlier mentioned spinning nature of the 
sound excitation at the inlet exit plane, the boundary condition describ- 

I 

ing the normal acoustic velocity, at this plane can be '-.•pressed in 
the following form 

f(r) = f(r) + if(r) (>) 

exit 

where the complex quantity f(r) describes the radial dependence of the 
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source. Using Eq. (5)» Eq. (9) can be rewritten as follows 


“ ?(r) 

(10) 

exit 


- f(r) 

(11) 

exit 



The sound excitation at the inlet exit plane could also be prescribed 
in terms of a spinning pressure wave expressed as follows 

^exlt “ 9^^^® ' 


where g(r) represents the radial variation of the pressure source. 
Suppressing the t and 9 dependence In Eq. (12) and separating into real 
and imaginary parts yields 


p«it ■ 


For the study of plane velocity wave propagation (i.e., m=0) the 

condition f(r) = a constant applies and similarly for a plane pressure 

wave propagation the condition g(r) = a constant applies. For a more 

general excitation, an appropriate combination of higher order spinning 

1 3 

modes, (i.e., Bessel functions ) must be used to describe the sound 
source. 

At a hard wall of the duct boundary, the appropriate boundary 
condition is 

<t>'^ = 0 (15) 

Equation (15) expresses the fact that fluid can not penetrate the wall. 
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To absorb the emitted turbomachinery noise in the front and 
aft ends of Jet engines* various types of sound absorbing liners have 
been developed by the aircraft industry. For example, point reacting 
liners which are essentially a series of sharply tuned Helmholtz 
resonators and bulk reacting liners which provide a broad band noise 
reduction by virtue of their porous lining have been developed. Though 
the point reacting liners are effective absorbers of sound over a 
narrow range of frequencies, they are better suited to sustain the 
operating conditions of a jet engine compared to the bulk reacting 
liners. To estimate the liner performance, the appropriate boundary 
conditions need to be applied at the interface of a point reacting 
liner. A good deal of controversy exists in the literature about the 
correct boundary condition to be applied across the liner in the pres- 
ence of a grazing flow, as to whether particle displacement or particle 

velocity normal to the liner is continuous. Matched asymptotic expan- 

4 15 

Sion studies by Eversman and Beckemeyer and Tester indicate that In 
the limit of zero boundary layer thickness, the results using the shear 
flow model approach the no shear model employing the condition of par- 
ticle displacement continuity. Nayfeh et al.^^ confirmed the conclusions 
of References 14 and 15 numerically. 

A physically meaningful explanation of the continuity of particle 
displacement is as follows. Consider an interface between two fluids 
being in relative motion. When a sound wave is incident on this inter- 
face from one of the mediums, it generates ripples on the interface and 
also creates a reflected wave In the original medium and a transmitted 


f 
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wave in the other median. The interface is a stream surface common to 
the flows in the two media. The boundary conditions at the interface are: 

(a) The pressure must be equal on both sides, and 

(b) The stream surface slopes must be equal on both sides, if S(r^,t) 
is the equation of the interface which deforms with respect to time, 
the statement corresponding to condition (b) is 


Dt 


“ 0 


(16) 


where is the substantial time derivative or time derivative fol- 

lowing the fluid particle. Equation (16) states that the adjacent fluid 
particles on either side move along the interface which indeed is the 
statement of particle displacement continuity. For the case of a rela- 
tive motion across the interface, continuity of particle displacement 
does not imply continuity of acoustic particle velocity, since acoustic 
particle velocity is the time derivative of particle displacement fol- 
lowing the fluid particle, it receives contributions from the convective 
terms which are different on either side of the interface. 

Hence the appropriate boundary conditions across the liner 
interface are: 


p' = Pp (17) 

C = C (18) 

P 


where C is the acoustic particle displacement and subscript p refers 
to the liner side. The equivalent condition to Eq. (18) in the continuim 


is Eq. (16) which states that any particle on the interface S(j^,t) 
remains on it for all times » i.e.. 
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11+ (V-V)S - 0 


and 


at S “ 0 


II + (V -V)S « 0 
3t 


(19) 


For the special case S(£, t) ■ n - C(?)e where n is outward normal 
coordinate and s is the coordinate along the boundary of the wall, Eqs. 
(19) become 


ii 

3t 

9S 


+ (rf> +4' )-~ + (4 ■*■4 )"^ (4n 4 ’o)’^7T ~ 0 

'r_ T.» >sc 't— t — 'tq 


S' 3S ""n ’’n'9n 

9S . /T . . * \ 3S . /7 . . • \ 9S 


( 20 ) 


+ (I +A.) ^ (A + 4 ,') ^+4 +^') -^=0 

3t ^s' ^s p 9s ^n ^n p 9n ^0 ^0 p 90 


Since the normal and azimuthal components of the mean flow are zero 
and the mean flow is assumed to be identically zero on the liner side, 
the above equations yield 


• \ "iwt - ^(s) -icut _ n 

10 ) 5 (s)e - 4g — e r - 0 


io)5 (s)e + 4 ^ = 0 

^ P 


( 21 ) 


The specific acoustic impedance of the liner is defined by 

I 

1 P, 


^£= 7 ^ 4 * 

c p 

0) 0) p 


(22) 
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where p and c respectively describe the values of density and sound 

u UJ 

speed at the wall. Substituting Eqs. (17)* (18), and (22) into Eqs. 
(21) yields 


+ I 


P c Z 

ui ui I 


Hi 


i_ (--E. 

Ss - 


P C z, 

u (U X 


(23) 


it is assumed that p c 2. is a piecewise constant function of s. Hence 

(I) to I 

Eq. (23) becomes 


P c Z„ 
w u £ n 


P' + 


ll l£l 

u) 8s 


(2i») 


As seen above, the condition of particle displacement continuity and 
that of particle velocity continuity are Identical If there Is no rela- 
tive motion between the two fluids on either side of the liner (i.e., 
ig 2 0). 

For a general liner Impedance p c Z = 0 + ix , Eq. (24) can be 

. Ui (0 £ £ X. 

rewritten by equating real and imaginary parts separately as 




^s ^ 

" P ■ 


(25) 


+ 6^ 
n ’^n 


P 


+ ii iE 

u) 3s 


(26) 


Decomposing velocities along and normal to the wall and noting that the 
normal component of the mean flow velocity is zero at the wall, Eq. 

(8) after separating into real and imaginary parts becomes 

p = -plo)^ + (27) 

p = ^ 5 ] 


(28) 
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Substituting £gs« (27) sikI (28) into Eg$> (25) snd (26) yields the 
following desired lined waii boundary condition in terms of the 
velocity potential and its derivatives only: 

i 

®^n‘^^n ” ♦s^s^ ^ ^ss^s"*” ^s^ss^ 


/2D3c^Sf~ --'i 

■!■ ( ■ ■*“ — [“ojA + 4 A 1 

Y“i ; 9s U) ^ ^s^s^ 


(29) 


and 


x4 +-6A 
^n ^n 


-p[-u)4 + A^Ag] - p 


— [(i)A 
0 ) s 


^ A + A A ] 
ss^s s^ss 




“3s (i) 


c 


[o)A + 


s’^s 


(30) 


Due to the complex nature of the reflection process at the 
inlet entrance plane, the precise form of the boundary condition at 
this location is currently not known. Rice^^ has argued that except 
for modes near cut off frequencies the assumption of no refiection of 
"internal" duct waves at the inlet entrance plane Is a reasonable one. 
As the primary objective of the present analysis is the development of 
the needed solution technique, the inlet entrance plane boundary condi- 
tion in the present study is specified in the following form 

P' =-p c 2^ Az (31) 

where p and c are the local steady flow density and sound speed at 
the entrance plane and is the local impedance. The analytical solu- 
tion for the propagation of a single acoustic mode, with cut-off fre- 
quency B, in a cylindrical or annular duct with constant mean flow Mach 
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number M is known. This solution provides the following Impedance 
condition for transmission without reflection: 


Z 

e 


U) * ~ 


(32) 


and values of S are available in Reference 13* For plane waves, the 
cut-off frequency, s, equals zero and Eq. (32) reduces to the more 
fami 1 iar resui t 


Since in the inlet case the steady flow velocity is not uniform at the 
inlet entrance plane, there will be a partial reflection of the principle 
mode. Furthermore, any additional modes excited by duct cross sectional 
area variations and steady flow velocity gradients will be also partially 
reflected. 

A general radially varying radiation impedance condition at the 
inlet entrance plane can be described by 

Z (r) = 9 (r) + ix (r) (33) 

e e e 

Application of Eqs. (8) and (33) to Eq. (31) yields the following, 
equivalent boundary conditions: 

<t>^c Xg + 0^(c 0g - ^^) = -0)1 + 


(35) 
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CHAPTER m 
HETHODS OF SOLUTION 

3.1 The Inlet Steady Flow 
3.1.1 Analysis of the Incompressible Potential Flow 

Analytical solutions for the nonlinear equations describing 
the steady compressible flow in ax i symmetric, axially nonuniform 
passages are not generally available and complex numerical solution 
approaches must be employed to obtain the desired flow description. 
Since the main objective of the current study is the description of 
the acoustic flow field, for which the steady flow is needed as an 
input, an approximate solution was used to obtain the needed steady 
inlet flow description. The approximate steady flow computation con- 
sists of a potential flow solution with a correction accounting for 
compressibility effects. An integral solution technique was used to 

compute the inlet potential flow utilizing a computer program developed 

18 1 9 

earlier at Georgia Tech, and Lieblein's correction ^ was utilized to 
account for compressibility effects. 

The incompressible potential flow is governed by Laplace's equa- 
tion and is subject to boundary conditions specifying the magnitude of 
the velocity normal to the inlet. The component of velocity normal to 
the solid surfaces of the inlet must be zero. A finite velocity distri 
but ion may be prescribed at the Inlet fan plane as the forcing boundary 
condition for this boundary value problem. The governing equations in 
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ii doubly connected region, R bounded by the inlet surface, B and a 
spherical surface of infinite radius, B' are 

- 0 In R (36) 

''♦■SIb - f - • S|b >37) 

and 

Iv^lg, - 0 (38) 

where 0 is the velocity potential created because of the inlet, n is 
unit normal vector pointing away from region R and F is the prescribed 
velocity normal to the inlet surface B. 

20 

The method developed by Smith et al. to solve the linear 
incompressible potential flow problem for bodies of arbitrary shape 
consists of transforming the differential equation (36) subject to the 
boundary conditions (37) and (38) into a linear surface integral equa- 
tion. In this method, the flow field created due to the inlet is 
thought of as due to a continuous distribution of sources and sinks on 
the inlet surface. Hence, the potential at a point p due to the source 
distribution on the inlet is 



where o(p‘) is the unknown source strength distribution on the inlet 
surface, B and r(p,p') is the distance between points p and p'. It 
can be easily verified that the velocity potential given by Eq. (39) 
satisfies Eqs. (36) and (38) for a finite value of a(p'). The source 
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distribution o(p') is dotermined by requiring it to satisfy Eq. (37) 
which leads to the formulation of the following surface integral equa-* 
tion'®*^’ for a 

2iro(p) - / o(p') Vl- ^(p*p, y] * n dB ■ F(p) - • n|p (40) 

The first term in Eq. (40) is the velocity normal to the inlet induced 
at p by the source at p'. The second term is the normal velocity compo- 
nent at a point p due to the source distribution over the remainder of 
the inlet surface. 

For flows past axisymmetric or two dimensional bodies Eq. (40) 
reduces to a line integral equation in a single plane. Furthermore, a 
numerical solution to the line integral equation is sought by assuming 
that the inlet surface is made up of many small straight line segments 
and that the source strength o is constant over each segment. The above 
discretization procedure leads to a set of linear simultaneous algebraic 
equations for source strengths over each segment which can be solved by 
standard matrix routines. The velocities on and off the inlet surface 
are calculated from the computed sou.-ce distribution (see Reference 18 
for a detailed mathematical development). 

3.1.2 Superposition of Solutions 

Since the Laplace equation (36) is linear, the principle of 
superposition can be utilized to obtain solutions for various free stream 
velocities and Inlet flow rates from the following two basic solutions 
(see Figure 2): 





Ik 


Let the computed velocities for Case 1 be denoted by (-q-’) and 

1 

U 

for Case 2 by {q-) , since V- is the reference velocity for Case 2. The 
F 2 

principle of superposition for these two cases yields 


f -Mf) 

CD I OB P 4 


(^1) 


The constants A and B are obtained by using the boundary values of y- 

m 

at Z ■ -« and Z ■ Zp. 

That is. 


u ^ C 

Urn ■^■l“Axl-*-B’ (-g~) X 0 ■ A 

« eo OB 


(^♦2) 


Since (tj-» tends to zero as Z-*-co as the velocity created at Z ■ 

f 2 

due to a suction at the fan Is zero. And at Zp, 

V V 

f I -1x1 + e(f) X I - ^ 


or 


B---V7 


(43) 


V F 

Since from Case I (y-) at Zp is Substituting Eqs. (42) and (43) 

CD 1 on 

into Eg. (4l) one obtains 


f ■ (f> * <r- 


(44) 


The ratio of the two reference velocities is determined by the pre- 

€C 

scribed free stream velocity and the mass flow rate through the fan. 
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3.1.3 Compress ibiUty Correction for Flow 

Next, since the inlet will be operating at high subsonic Mach 
numbers during takeoff or landing configurations, the incompressible 
solution obtained earlier needs to be corrected to account for compres- 
sibility effects using the semi- empi r leal equation proposed by Lieblin 
and Stockman. The justification for such an ad Njc approach to arrive 
at a compressible mean flow description within the inlet is that obtain- 
ing solutions to the exact nonlinear compressible potential flow equa- 
tion in a complex duct shape is beyond the scope of the present study. 

The empirical equation developed In Reference 19 to obtain "com- 
pressible" velocity, from incompressible velocity, V. at a point Is 

V./V. 

P. I I 

V • V. (-i-) (^5) 

c . i 

^c 

where V. is the area weighted average incompressible velocity across the 
duct at a given station, pj is the incompressible density which is also 
equal to the stagnation density and is the average compressible den- 
sity across the flow passage. If the duct were a circular cylinder, the 
exponent Vj/V. equals unity and Eq. {kS) is a statement of continuity. 
However, because of area variations in an aircraft inlet Vj/V. does not 
equal unity and it expresses the Influence of geometry on the "compres- 
sible" velocity distribution. 

The unknowns and in Eq. (^5) are related through the isen- 
tropic gas dynamic relations as follows. Principle of mass conserva- 
tion at an axial station requires that 
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m) 


Dividing Eq. (i»6) by the critical velocity V” and using the gas dynanic 
(23) 

relations one obtains 


V, p, Y-l 1/2 

.1 D. P: 


i^7) 


V. 

The ratio — x can be related to the free stream Hach number by making 
V" 

use of the ^sentropic gas dynamic relations again and one can show that 


V. 

I 



c 



OD 
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/ I + 


Y-1 
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00 


(^8) 


Vi 

Denoting for simplicity — as a and — as B equation (1»7) reduces to 


the following nonlinear algebraic equation for B: 
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(— ) 
4+r 


2 

a 


0 


(^♦9) 


Equation ikS) is solved by the classical Newton- Ra 1 phson scheme as 
fol lows: 

Setting 

e/ ^ 2 Y+1 /Y-K 2 

P(b) = 8 - 8 - 

one obtains 

F* = ■— - 28 - (y+Db^ 


The iteration scheme is 


F(& ) 

- PTT’ <5<" 

n 

A reasonable guess for the initial value, 8^ is 0.9 and the iteration 
scheme converges very rapidly to a physically meaningful value. Know** 
ing $ the conpressible velocity is calculated by using Eq. (A5) 

3.2 Finite Element Solution of the Acoustic Equations 
Due to its apparent advantages, the application of the FEM in 
the solution of a variety of engineering problems has been rapidly grow- 
ing in recent years. A detailed discussion of the FEM can be found in 
Reference II. Due to its suitability for handling problems Involving 
complex geometries and mixed boundary conditions the FEM is used in the 
present investigation in the solution of the Inlet wave equations. 

Any numerical procedure developed to solve a problem in the realm 
of continuum mechanics essentially converts the problem with an Infinite 
number of degrees of freedom to one with a finite number of degrees of 
freedom. The solutions so obtained are expected to approach the exact 
solutions as the number of degrees of freedom is increased. In the 
finite element method the continuum region under consideration is sub- 
divided into a number of elements and the variation of the flela varia- 
ble within each element is prescribed in terms of its value at a number 
of preassigned points in each element (I.e., nodes) and a set of known 
Interpolating functions. Once the values of the field variables at the 
nodes are determined, the behavior of the field variables within each 
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el^nent and hence in the overall continuum is known. Thus, the probl^ 
of determining the field variables at an infinite nunber of points is 
reduced to one of finding their values at a finite number of nodes. 

To obtain the needed solutions, Eqs. (6) and (7) which describe 
the wave propagation in ths inlet have been transformed into integral 
equations utilizing the Galerkin Method. The resulting integral 
equations were then solved using the FEM. The solution involves the 
following six operations: 

3.2.1 Subdivision of the solution domain: In the finite element method 

the commonly employed element shapes for discretizing the domain of 
interest are triangles, rectangles or i ioparan^tric triangles (i.e., 
triangles with curved edges) for two dimensional or axisymmetric problems 
and tetrahedrons, cuboids or isoparametric tetrahedrons for three dimen- 
sional problems. The choice of the element shapes in solving a particu- 
lar problem is determined by the geometry of the domain of interest and 
the level of complexity that can be employed. For two dimensional prob- 
lems, straight edged triangular elements are by far the most cotmionly 
used elements because of the relative ease in closely approximating com- 
plex shaped domains as compared to rectangular elements and the lower 
level of complexity in the finite element methodology as compared to the 
isoparametric triangles. Hence, straight edged triangular elements 
with three or six nodes have been chosen in the present study for sub- 
dividing the inlet duct. A computer code has been developed which sub- 
divides a duct into triangles with three nodes and catalogues the 
geometric nodal locations, nodal numbers and element numbers. The assent»ly 
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of nodes and trianyles is such that each element number corresponds to 
three nodal numbers and each nodai number is associated with a fixed 
number of elements whose number may vary between two and six depending 
upon the node location (see Figure 3)* For the six node elements, the 
vertices and the midside points of the triangle are normally chosen as 
the nodes. in such a case, the node and element numbers are catalogued 
so that each element number corresponds to six nodal numbers and each 
nodal number is associated wi h a fixed nunber of eluents varying 
between one and six (see F<i:ture 4). 

3.2.2 Selection of Interpol ht ion Functions: Polynomials are the popu- 

lar interpolating functions because of their desirable mathematical 
properties of completeness and ease in mathematical manipulation. The 
variation of the field variable (acoustic potential in this case) 
within a three node triangle can be expressed in terms of a linear 
Lagrangian polynomial in two dimensions. Hence, three node triangles 
are also known as linear triangles. A quadratic representation of the 
acoustic potential can be obtained using six node triangles since a 
quadratic Lagrangian polynomial in two dimensions has six coefficients. 

A linear transformation from the global Cartesian system to a 
local "area coordinate system" for each element has been found to simplify 
the mathematical operations involving triangular elements considerably. 

In this transformation, a triad set of coordinates for each element, 

(e) (e) (e) 

(L. , L. , L' ') out of which only two are linearly independent 

I J K 

replaces the diad set of coordinates, (r,z) through the following rela- 


t ionsh i ps : 
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l' Z, * Lj* Z. * L* Z, 


. e . 1 e , e 
Lf r. + Tj + 


L? + L? + L* 
I j k 


( 51 ) 


Solving Eqs. (50 for L®, L? and L^, 


U ■ (aj + bjZ c.r)/2A 


Lj - (aj + bjZ + Cjr)/2A 


(52) 


where 


'•k ■ * '=k''>^^* 


j det 


I Z. r. 


. Z. r. 


' ^ % 


Area of triangle 


and 


a. ® z.r, - z, r.: a. « z, r. - z.r a, = z.r. - z, r- , 
I jk kj j ki ik k ij j' 


’ '•j ■ 

b. = r, - r.; 

J k 1 

b, = r . - r . 
k 1 J 

- ^ - ^j-> 

c. = 2, *• Z. ; 
J i k 

c. “ Z. - Z. 
k J 1 


For a triangle with three nodes at the vertices, the interpolating 


functions are simply the area coordinates. Thus the acoustic potential 
over the e^^ element varies linearly and Is given by 


M 
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^(r,z) - L*(r,2)0j Lj(r,z)^j + L®(r,z)|^ (53) 

i(r,z) - L*(r,z)$^ + Lj{r,z)Jj + L®(r,z)0j^ (5^) 

A quadratic variation for the acoustic potential over the element Is 
achieved by choosing triangles with six nodes, three being corner nodes 
and the remaining three being mId-sIde nodes (see Figure 5) and is 
described by the foIlcM/Ing relationships: 

5^(r,z) « N®(r,z)^j + N?(r,z)^j + N®(r,z)|j^ + N®(r,z)^^ 

+ N®(r,z)| + N®(r,z)J (55) 

m m n ’^n 

$(r,z) - N®(r,z)J. + N?(r,z)^j + N®(r,z)$j^ + N®(r,z)^^ 

+ N*(r,z)i + N®(r,z)* (56) 

m n ^ n 


where the quadratic Interpolating functions are related to the area coor 
dinates L®, lJ, and L® by the following: 

For Corner nodes. 


N® = 2(Lp^ 

nJ - 2(L?)^ 

N ® = 2 ( l ®)2 
k k 


L? 

I 


L? 

J 


For mid-side nodes 


N=-i.L?L? 
Z I J 


m 


l*L?L® 
J k 


i»L®L? 
k I 


(57) 


( 58 ) 
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3<2.3 Establishing Etenental Relations: The Galerkin method is 

applied to the governing differential equations and boundary condi- 
tions to develop the matrix equations that express the properties of 
individual elements in terms of the unknown nodal values. Denoting 
either of the partial differential equations, 6 and 7f by the operator 
application of the Galerkin and FEM methods^ ^ yields the fol- 
lowing relationships: 

l III LU')dvl'* - 0 m-1,2 N (59) 


where the integration is performed over each element. In Eq. (59) N 
is the total number of nodes in the problem under consideration and E 

is the number of elements, it should be noted that N is zero for all 

m 

elements not having the nodal point, m, as a vertex. Equation (59) pro 
vides 2N equations for the 2N unknown nodal values. 

In evaluating Eq. (59) over each element the following integral 
involving the area coordinates and arbitrary exponents a,b, and c is 
needed: 


p a A b 


// (Lp (Lp (Lp drdz 


albic! » (e) 
{a+b+c+2)! 


(60) 


(see Appendix A for the proof of Eq. (60)). 

There exists^ ^ a mathematical restriction In the choice of 
interpolating functions stating that the values of <j>* and its partial 
derivatives up to the highest order appearing in Eq. (59) must have 
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representat ion as the element size shrinks to zero, if linear inter- 
polating functions are chosen, the second order partial derivatives of 
will be identically zero. To avoid this, all terms in Eq. (59) con- 
taining second order partial derivatives of (p* can be reduced to terms 
containing first order derivatives of using Green's theory for a 
plane geometry (also known as integration by parts). The boundary 
conditions, given In the previous chapter, are introduced into the 
boundary integrals that are obtained as a result of the above-mentioned 
integration by parts. However, for an Inlet lined with point-reacting 
liners carrying a steady mean flow the boundary conditions at the 
lined wall, as given by Eqs. (29) and (30), involve second derivatives 
of . Hence, linear interpolating functions cannot be used since the 
terms containing second derivatives of ^ will be identically equal to 
zero. In such a case, one has to choose quadratic or higher order rep- 
resentation of 0' over each triangle. In the present study, quadrat •; 
interpolating functions have been used to handle the lined wall boundary 
conditions. A point to note Is that If there is no mean flow in the 
duct, the lined wall boundary conditions do not contain any second order 
derivatives of p* and again one does not need a quadratic representation 
for <p to treat the no mean flow case. After the choice of interpolating 
functions is made, five different elemental relations are developed 
depending upon the location of the triangles under consideration; that 
is: (I) triangles interior to the flow region; (2) triangles adjacent 

to a hard wall; (3) triangles adjacent to a lined wall; (4) triangles 
adjacent to the fan plane; and (5) triangles adjacent to the entrance 


38 


plane (see Appendix B for the mathematical development). 

3 . 2 .I 1 Assembly of Element Equations to Obtain a System of Algebraic 
Equations: Based on the ordering system defined in the first step, 

the individual el«nent equations are ccMnbined Into a matrix equation 
describing the properties of the potential in the domain under 

study. Since each node is only affected by adjacent elements, the 
resulting matrix is banded. 

3 . 2.5 Solution of the System of Equations: Considerable amount of 

information Is currently available about the solution of large, 

banded matrices. in tiw present study, CDC subroutine 6LSWNP is used 
25 

for calculations. The BLSWNP subroutine has been developed to solve 
efficiently a system of linear algebraic equations with a banded coeff! 
cient matrix. Let kl and k2 be the bandwidth of the lower and upper 
triangles of the NxN banded matrix, excluding the main diagonal. All 
the elements outside this bandwidth are equal to zero. In the BLSWNP 
Subroutine only the nonzero elements (i.e. at most N x (Kl+k2+l)) are 
stored such that the columns of the condensed matrix are the diagonals 
of the banded matrix and rows are stored as rows of the banded matrix 
(see Figure 6). Such a condensed storage scheme reduces the computer 
memory requirements. Further, the solution of the system of equations 
is obtained in two sequential steps as described below. 

To solve the banded system of equations, 

Ik] {X} - (b) (61) 


decompose the problem into 



ILJ {y} - {b} 


i<0 

(62) 

and 

[U] {x} - {y} (63) 

where [L] and [U] are lower and upper triangular matrices derived 
from [K] and {y} is an intermediate solution vector. Since solution 
of triangular matrices can be obtained without pivoting, the above 
scheme is computationally quite efficient. 

3.2.6 Additional Calculations: Once the velocity potential is obtained 

at the nodes, additional variables such as the axial and radial 
velocities and acoustic pressures can be calculated at any point in 
the domain, by utilizing the computed solutions of 41 ' together with 
Eqs. (53) and (54) for a three node element case or with Eqs. (55) 
and (56) for a six node eluent as the case may be. 

The results obtained can be checked for the recovery of the 
boundary conditions and conservation of acoustic energy. The time 
averaged acoustic intensity in the a^^ direction for an irrotat ional , 
uni form- entropy flow to second order in acoustic quantities is given 

26 

by the following expression (also see Appendix C for the derivation) 

V(j) ‘ a 2 (Viji ■ a ) V(|» 

<j’> . a = <p'V;|,'> • a + _ <p' > + — .<p'v 4 >'> 

PC c 

+ P ( 74 ) •<74i'74>'^)'a (64) 

where J is the acoustic intensity vector, a is unit vector in the 
direction along which acoustic intensity is needed, and < > describes 
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time average over a tong period of time. The time average energy flow 
across a surface S is then 

<E> ■ / <1 ' > • n d S {65) 

S 

where n is the unit outward normal to surface S. 

in the present study, the time average acoustic energy flows at 
the entrance to the duct, at the exit plane of the duct and at the duct 
wails have been computed using Eqs. (64) and (65) (see Appendix C for 
a detailed finite el&i^nt evaluation procedure). The effectiveness of 
a liner in absorbing the sound is expressed by the decibel reduction 
defined by 

dB . =lOlog,rt 

reduction ^10 

As a check on th«* arniracy of the developed solution, the difference 
between the energy flow into the duct and out of the duct should equal 
the energy absorbed by the lined walls; i.e., acoustic energy should be 


<E>. 


I nput 


<E> 


output 


( 66 ) 


conserved. 


To check the accuracy of the developed FEH computer progran, 


solutions for the problems of plane and spinning wave propagation 
through a hard walled annular duct with constant mean flow have been 
obtained for comparison with available analytical solutions. For this 
case, Eqs. (6) and (7) reduce to the following form: 


r 

— * ‘ 4> + 2Mu>«^z = 0 (68) 

r 

where w Is thw frequency (non-d imensional I zed by the outer annulus 
diameter and the sound speed), m is the spinning mode number and M Is 
the constant mean flow Mach number. 

The hard wall boundary conditions are described by the following 
expressions: 


= 0 
“ 0 



0 (Inner wall) 

1 (outer wall) 


(69) 


and the sound source boundary condition at the exit plane (i.e., z= L) 
is given by 


Assuming no reflection at the duct entrance plane, where z » 0, the 
following boundary condition applies: 

p “ ; p = -2^2 (71) 

The needed expressions for the impedance Z and f(r) are given below. 

The exact solution to this problem is given by 

i sin[k(z-L)] (72) 

i cos[k(z-L)l (73) 

For the plane wave case (i.e. , m > 0), the various quantities in 
Eqs. (70) through (73) take on the following form: 

f(r) = 1 ; Z = 1 ; k = (74) 

The real and imaginary components of the velocity potential as calcu- 
lated by the FEM for an annular cylinder with a = 0.5 and L = I are 
compared with the exact values computed using Eqs. (72) and (73) in 
Figure 7 for the case of M « 0.5 and w = 2.0. A good agreement between 
the FEM and analytical solutions is shown; good agreement has also 
been obtained when the predictions for the acoustic velocities and 
pressures were compared. 

The expressions describing the propagation of a spinning wave 
with lobe number m ar.d radial mode p are: 






f(r) 


C ten 

my 


<D - B^O - M^) 

u)M + /w^ - B^(l - M^) 


(75) 


and 

Mh) ■<■ i/u)^ - (l ~ M^) 

I - 

where and 8 in Reference 13) are tabulated in Reference 13* 

my mu 

Results of FEM calculations for the velocity potential in the 
previously described annular duct with m=4, y=0, w=6 and M = .5 
are presented in Figure 8. The analytical solutions, given by Eqs. 

(72) and (73) are also shown in Figure 8 and good agreement between the 
two solutions is noted. Similar comparisons for other acoustic varia- 
bles such as acoustic velocity and pressure have also shown good agree- 
ment . 

For the FEM calculations shown above, the duct was subdivided 
into 220 3“node triangles with 136 nodes. For a Mach nunber of 0.5 
good agreement between FEM calculations and the exact solutions was 
obtained for values of w up to 10. At this point, (u =10, M = 0.5) 
there are about 13 nodes per wave length parallel to the annular cylin- 
der axis. Clearly, the accuracy lost at higher frequencies can be 
recovered by a finer element subdivision or a more elaborate description 
of the dependent variable <p within each element (e.g., quadratic repre- 
sentation of d) in each triangle). 
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4.2 Comparison of Acoustic Calculations by Using 
Linear and Quadratic Triangular Elements 

The problem of plane wave propagation in the annular cylinder 
carrying a uniform steady flow considered above is chosen to cwpare the 
accuracy obtained by using linear and quadratic elements for various 
frequencies. For the case of linear elements, the cylinder is divided 
into 140 3"node triangles with the total number of nodes being 94 
(see Figure 9) while for the case of quadratic elements, it is divided 
Into li»0 6-node triangles with a total of 327 nodes (see Figure 10), 94 
out of which are the corner nodes and the remaining 233 nodes are the 
mid- side nodes. 

FEM calculations have been performed for 5 angular frequencies 
(viz., u) * 1 , 5, 10, 15, and 20) using linear and quadratic elements. 

The amplitude and phase of the axial acoustic velocity obtained by FEM 
programs are shown in Figures 11 through 20 along with the exact solu- 
tions derived from Eqs. (72) and (73)- For the case of w = i and o) = 5, 
the results obtained by both the FEM programs agree with the exact solu- 
tion to within 1 to 2%. For the case of w = 10, quadratic elements give 
results accurate to 4^ while as linear elements give results accurate 
up to S%. The predictions by linear elements for u) = 15 and 20 are 
about 30^ off the exact solution and they do not even show the proper 
trend. However, quadratic elements yield results which are at worst 7^ 
off for li) = 15 and 12^ off the exact solutions for w = 20. Moreover, 
even at n) * 20 the results of quadratic elements Indicate the proper 
trend (e.g., the acoustic axial velocity amplitude remains reasonably 
constant). These results unequivocally prove that quadratic representation 
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Comparison of Axial Acoustic Velocity Amplitude 
for Plane Wave Propagation In an Annular Cylinder 
as Predicted by Linear and Quadratic Elements 
(u) ■ 1.0, Mach « 0. 5) 
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Figure 16. Comparison of Axial Acoustic Velocity Phase for Plane 

Wave Propagation in an Annular Cylinder as Predicted by 
Linear and Quadratic Elements for co = 10. 0 and Mach = 0. 5 
(Symbols defined in Figure 11). 




Figure 18. Comparison of Axial Acoustic Velocity Phase for Plane 
Wave Propagation in an Annular Cylinder as Predicted 
by Linear and Quadratic Elements for w = 15,0 and 
Mach = 0.5 (Symbols defined in. Figure 11). 
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Figure 13. Comparison of Axial Acoustic Velocity Amplitude for Plane 
Wave Propagation in an Annular Cylinder as Predicted by 
Linear and Quadratic Elements for w = 20.0 and Mach = 0.5 
(Symbols defined in Figure 11). 
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of ^ over each element yields much better results compared to linear 
representation for higher frequencies, for a fixed number of elements. 


A. 3 Comparison Studies for Lined Wall Annular Cylinders 
In order to evaluate the accuracy of the finite element method 
and to de-bug the computer codes, the case of sound propagation in an 
annular cylinder with uniform steady flow and lined walls has been inves- 
t i gated. 

Because a comprehensive set of results including radial and 

axial profiles of acoustic variables along with the dB . for the 

reduct i on 

case of a cylindrical duct is not available, comparison was made with 

2 

the results of Baumeister for the case of a rectangular duct carrying 

a uniform steady flow. A rectangular duct is a good approximation to 

the annular cylinder if the radius of curvature of the cylinder is large 

and if the radius ratio is sufficiently close to 1. The investigated 

problem is example #3 in Appendix F of Reference 2. The geometry and 

associated parameters are shown in Figure 21. The radius of the annular 

cylinder was increased from 100 to 1000 to check the convergence of the 

dB . to a constant value. The dB . obtained is 4.726 which is close 
red red 

to the 5.6 value obtained by Baumeister. Baumeister has solved the 

28 

same problem by using the generalized wave envelope transformation. 

A new variable, ijf, defined as 
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<J^(r, 2 ) * e (r,z) 


(76) 


where 6 is a free constant, is substituted into Eqs. (67), (68), (70) 
and ( 71 ) to obtain a new set of partial differential equations and 


— Plane Pressure, p 



0 z 


Figure 21. Annular Cylinder Geometry Siniulating Rectangular Duct. 
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boundary conditions for the variable ij;. Having chosen 6 to be equal to 
the appropriate wave number, Baumeister has obtained a dB^^^ of k.072. 
The obtained by the FEM program lies in between the two values 

obtained in Reference 2. Figures 22 through 25 indicate the variation 
of acoustic pressure (I.e. , magnitude and phase) as predicted by the 
FEM program and Reference 2. Also plotted are the corresponding hard 
wall duct solutions for the same flow conditions, for the purpose of 
indicating the effect of the liner in reducing the pressure amplitude. 
With the exception of the results in Figure 22, the agreement between 
the two sets of computations is good. Since the liner is known to 
attenuate the sound wave, one would expect a reduction in acoustic 
pressure magnitude as one moves from the fan plane towards the entrance 
plane of the duct. While the FEM results agree with this intuitive 
argument, the results of Reference 2 indicate an opposite trend. A 
possible explanation for this is as follows. The finite difference 
grid used in Reference 2 employs 100 points to model the region of 
interest while as the finite element model employed 32? points to model 
the same region. Hence, one could anticipate a better agreement between 
the two solution schemes if the finite difference grid is made finer. 
Figure 26 shows the variation of acoustic power absorbed by the lined 
upper wall. The total acoustic energy absorbed by the liner should 
equal the difference between the acoustic energy at the entrance and 
exit planes. This acoustic energy balance has been checked and acoustic 
energy is conserved with an error of 8^. The prescribed radiation impe- 
dance condition at the exit plane and the pressure boundary condition 



0 0.25 Z 0.5 


Figure 22. Comparison of Acoustic Pressure Magnitude Along 
the Upper Wall of the Rectangular Duct as 
Predicted by FEM Program and Reference 2. 

Hard Wall Solutions Shown in Circles. 






Figure 2k. Comparison of Acoustic Pressure Amplitude 
Along the Lovver Wall as Predicted by FEH 
Program and Reference 2. Hard wall solution 
shown in circles. (Symbols defined in 
Figure 22.) 
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Figure 25- Comparison of Acoustic Pressure Phase Along 
the Lower Wall as Predicted by FEM Program 
and Reference 2. Hard wall Solutions 
shown in circles. (Symbols defined in 
F Igure 22. ) 
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at the entrance plane are recovered with a maximw error of 6^. Errors 
obtained in other computed soft wall cases are smaller than the above- 
n^ntioned 6 rnd 8 percent. Computations for hard walled annular ducts 
produce acoustic energy conservation with an error of 10 and the 
radiation impedance condition and the pressure boundary condition are 

-ii 

recovered with 10 % error. 

4. it Comparison Studies wi ^ the Integral Technique 
for the QCSEE Inlet 

The problem of prescribing the correct radiation impedance con- 
dition at the entrance plane of an inlet carrying a variable mean flow 
has not yet been solved. However, for the case of no mean flow, Bell, 

3 

Meyer and Zinn have solved the problem of wave propagation in an inlet 
using the Green's function approach in the interior of the inlet and 
also in the far field. With the knowledge of the solution in the inter- 
ior and the exterior regions of the inlet a radially varying radiation 
impedance can be calculated at the entrance plane of the inlet. The 

calculated Z was incorporated in the FEM program that was specialized 
e 

to handle the no mean flow case and the solutions obtained by the FEM 
program and the integral approach for various frequencies and for both 
hard and soft walled QCSEE inlet one compared in this section. 

For the case of no mean flow {i.e., =0) Eq. (8) reduces 

to 


p‘ = ioj(|>' 

since, p = I as the density of air In the inlet equals the ambient air 
density. Thus, the acoustic pressure is directly proportional to the 
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acoustic velocity potential and comparisons have been made with the 
acoustic potential rather than the acoustic pressure. 

As described in Section 3*2.3> in the absence of mean flow one 

does not need a quadratic representation for the acoustic potential ^ 

to handle the lined wall boundary condition. Hence, the comparison 

studies described in this section have been performed using the 3~node 

tr iangul izat ion scheme of the QCSEE inlet (see Figure 3). in both the 

FEM and integral approaches the sound excitation at the fan plane is 

prescribed as a plane velocity wave of unit amplitude and a phase angle 

of l80°. Table 1 shows the distribution of Z at the nodes located at 

e 

the entrance plane of the hard walled QCSEE inlet as calculated by the 
integral approach for the plane wave propagation of angular frequencies 
(i) *= 1,2,5 and 10. The above radiation impedances were incorporated In 
the FEM program and the acoustic velocity potentials obtained by both 
the FEM and the Integral approach are plotted in Figures 27, 28, 29 and 
30 for plane wave propagation of angular frequencies to = 1,2,5 and 10 
respectively. One observes that the agreement between the two results 
is excellent and particularly so at the low frequencies of u> = 1 and 2. 
Though the results obtained by both methods are slightly different for 
0 ) = 5 and 10, they indicate similar trends in their behavior. Th.e inte- 
gral approach used only 50 source points on the inlet surface whereas 
the FEM used 212 nodal points to map the domain of interest. A better 
agieement between the two methoas at higher frequencies can be obtained 
if more source points are utilized In the integral approach to represent 
the inlet surface. 
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Figure 23 . Acoustic Velocity Potential Distribution Along the Inlet Upper Wall 
(u) = 5-0, No Flow, Hard Wall, Plane Wave Excitation; Symbols 
defined In Flcjure 27). 



Comparison studies have also been conducted for the case of plane 
wave propagation in the QCSEE inlet with lined wails. In the absence 
of mean flow the condition of particle displacement continuity and that 
of particle velocity continuity are identical. The specific acoustic 
impedance of the inlet upper wall is chosen to be = 0.l6-f0.3^, a 
representative value chosen from Reference 2 for the case of plane wave 
propagation in a lined rectangular duct carrying no mean flow. The 
center body is still prescribed to be a hard wall. The sound excitation 
at the fan plane is again prescribed as a plane velocity wave of unit 
amplitude and a phase angle of l80”. 

Table I shows the distribution of 2 at the nodes located at the 

e 

entrance plane of the soft walled QCSEE inlet as calculated by the 

integral approach for the plane wave propagation of angular frequencies 

oj * I and 2. The acoustic velocity potential distribution along the 

lined upper wall of the inlet obtained by the FEM and the integral approach 

for angular frequencies cj = 1 and 2 are shown in Figures 31 and 32, 

respectively. One observes that the agreement between the two methods 

for the lined wall case is also excellent. The dB . .. obtained for 

reduct i on 

both the frequencies by the two methods agree to four significant places 
\for u) = 1 dB . equals 24.879^ and for w = 2 dB . . . equals 

13 . 35^^)* Figures 33 and 3^ show the variation of the acoustic energy 
absorbed by the lined wall for u = 1 and 2 respectively. With this 
information the conservation of acoust’c energy is checked by comparing 
the inflow and outflow of acoustic energy in the domain of interest. 

It has been found that the acoustic energy is conserved to within 3^- 


♦rro! -91 ' 
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Acoustic Velocity Potential Distribution Along the Inlet Upper Wall 
(w = 1.0, No Flow, Soft Wall, Plane Wave Excitation, Symbols 
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Table 1. "Exact" Radiation Impedances from the Integral Technique 
for the QCSEE Inlet With No Mean Flow 


Z ■ 2^ iX for the QCSEE Inlet 
e e e 

Hard Wall Soft Wall 


Ui U) 

Node ~ ~ 

Number 1.0 2.0 5.0 10.0 1.0 2.0 


0.2629- 

0.4415- 

1.0942- 

0.8815- 

1.1780- 

1 . 7210 - 

i 0.5333 

i 0.4676 

10.3231 

10. 1020 

11.7072 

12.2895 

0. 2^*27- 

0.7204- 

1.1589- 

0.9365- 

0 . 7214 - 

1 . 5664+ 

\0.50Bk 

10.5629 

10. 1141 

1 0. 0672 

1 1 .2168 

10.7554 

0.26A2- 

0.9994- 

1.2237+ 

0.9914- 

0.2648- 

1.4119- 

i 0.6385 

10.6583 

1 0.0949 

10.0323 

10.7265 

10.7788 

0 . 3100 - 

1.2422- 

1.1573+ 

1.0679- 

0.2161- 

1 . 0600- 

i 0.7535 

10.6592 

10.0703 

1 0.0600 

10.6513 

10.9674 

0.3463- 

1.4175- 

1 . 1 046- 

1 . 0698- 

0.1965- 

0 . 9196 - 

i 0.8390 

i 0.6405 

1 0.0486 

1 0.0788 

1 0.o26') 

1 1.0234 

0.3724- 

1.5344- 

1 . 0754- 

1 . 0287 + 

0 . 1869 - 

0 . 6535 - 

i 0.8993 

i0.6l95 

1 0.1659 

1 0.0386 

10.6173 

1 1.0511 

0. 3880- 

1.6012- 

1.0625- 

1 . 1025 + 

0 1822 - 

0 . 8222 - 

i 0,9348 

i 0. 6044 

1 0.2440 

1 0.1934 

i 0.6136 

1 1.0655 

0.3932- 

1.6231- 

1 . 0588- 

1 . 1572 + 

0. 1808- 

0 . 8131 - 

i 0.9466 

i 0.5991 

1 0.2708 

1 0.2360 

i 0.6126 

1 1.0701 
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error for both the frequencies. 

Since the prescription of the "exact" radiation Impedance con- 
dition at the entrance plane of an Inlet carrying a variable mean flow 
is not possible with the present day knowledge, the simple "no reflection" 
radiation impedance condition, p c, was used in most of the investi- 
gated cases. The p c condition yields no reflection of the incident 
wave only if the propagation of plane waves Is considered and if the 
reflections from the inlet walls are not significant. 

in the absence of mean flow, if the inlet could be app* "ximated 

as a tube of constant cross sectional area, the problem of pijne wave 

propagation within such an inlet of finite length and in the far field 

2g 

can be solved by the classical Weiner Hopf Integral technique for 
frequencies below the cut off frequency corresponding to the lowest 
order spinning mode, namely the first radial nK»de. For the inlet 
modelled as a circular pipe the cut off frequency for the lowest order 
spinning mode is ■ 3.832. ^ Hence, for frequencies greater than 
plane wave solutions by the Weiner Hopf technique do not exist. The 
integral solution obtained by the Weiner Hopf technique yie’Hs a con- 
stant (in general complex) impedance condition at the entrance plane of 
the inlet to be henceforth denoted as "Weiner Hopf radiation impedance 
condition," whose value for w ■ 1.0 and 2.0 equals 0.2335 - '0.5555 
and 0.8317 ■ /O.6528 respectively. One could prescribe this radiation 
condition for the corresponding frequencies in place of p c condition 
for the case of no mean flow. 

Studies have been conducted to evaluate the Influence of the 
"exact" radiation condition, p c "no reflection" radiation condition and 
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the Weiner Hopf radiation impedance condition (where appiicabie) on 
the wave structure within the QCSEE inlet for the case of piane wave 
propagation of angular frequencies u ■ 1,2,5 and 10. Figures 35 and 36 
show the Influence of the above three radiation impedances on the 
acoustic velocity potential distribution along the hard upper wall of 
the inlet for « - I and 2 respectively. One observes In Figures 35 
and 36 that the wave profile predicted by prescribing the Weiner Hopf 
radiation condition matches quite well with the one predicted by pre- 
scribing the "exact" radiation condition. Thus at low frequencies the 
Weiner Hopf radiation condition can be viewed as an "average" constant 
value of the exact radiation condition across the inlet entrance plane. 

The wave structure obtained by imposing the "no reflection radia- 
tion condition is markedly different from the one corresponding to the 
"exact" radiation condition for w ■ 1.0 (see Figure 35) • However for 
the case of u ■ 2.0, the profile corresponding to the "no reflection" 
radiation condition tends to approach the wave profile corresponding 
to the "exact" radiation condition (see Figure 36). Figure 37 shows 
that the wave profiles obtained by prescribing the p c radiation condi- 
tion >f.id the "exact" radiation condition for a plane wave propagation 
of u) “ 5 metili quite well along the length of the inlet. An almost 
exact match between the two profiles is observed for the case of (j “ 10 
as shown In Figure 38. A physical explanation for this observation is 
as follows. The propagation of high frequency plane waves in a duct 
could be appfoximated a' ^he wave propagation in a ray tube where the 
reflected component is negligible, in which case, the "no reflection" 



igure 35* Influence of Different Radiation Impedances on the Acoustic Velocity Potential 
Distribution Along the Inlet Upper Wall (ui « 1.0, No Flow, Hard Wall, Plane 
Wave Excitation). 






igure 3b- Influence of Different Radiation Impedances on the Acoustic Velocity 
Potential Distribution Along the Inlet Upper Wall (oi * 2 . 0 , f4o Flow, 
Hard Wall, Plane Wave Excitation; Symbols defined In Figure 35)- 
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radiation condition indeed approaches the "exact" radiation condition. 
This fact is also evident from the values of the "exact" radiation irnpe* 
dances given in Table I for o) ■ 1,2,5 and iO. As u increases, one 
observes from this table that the "exact" radiation condition approaches 
the value of 1 + 01, which is the "no reflection" p c radiation condi- 
tion in terms of the nondimensional quantities. 

k.5 "Compressible" Mean Flow Calculations for the 
QCSEE and Bellmouth Turbofan Inlets 

To predict the acoustic properties of practical inlet configura- 
tions, the QCSEE (Quiet, Clean, Short-haul Experimental Engine) inlet^ 
and the Bellmouth inlet^^ have been chosen for the present study (see 
Figures la and lb). 

Before proceeding with the inlet acoustic analysis the b'^jhavior 
of the inlet mean flow must be determined. As explained earlier, the 
steady flow in the inlets is approximated using a solution for an incom- 
pressible steady flow with a compressibility correction. 

The incompressible solution was obtained by solving an integral 
formulation of Laplace's equation with the inlet boundary divided into 
a number of straight line segments. For the QCSEE inlet, the number of 
straight line segments equals HO and It equals 90 only for the Bellmouth 
inlet, because of the simpler geometry of the Bellmouth inlet compared 
to the QCSEE inlet. The chosen free stream and fan plane velocities 
correspond to a free stream Mach number of 0.12 and an average exit 

9 

Mach number of 0.52 for the QCSEE inlet and correspond to a free stream 
Mach number of 0.0 and an average exit Mach number of 0.52 for the 



Bet In»uth inlet. The resulting solutions of the Laplace's equation are 
superposed as described in Section 3.1.2 to yield the incompressible 
mean flow descriptions for both the Inlets. Application of the com- 
pressibility correction to the incompressible mean flow velocities as 
described in Section 3.1.3 yields the "compressible" velocity descrip- 
tion of the mean flow in the interior of the inlets. 

The mean flow velocity computations needed for the acoustic equa- 
tions have been performed at the nodal locations of the quadratic fin- 
ite element triangul izat ion schemes of the QCSEE and the Bellmouth 
inlets (see Figures ka and Ab), Tables D. I and D.2 in Appendix D con- 
tain the mean flow data at the corner nodal points along with the 
coordinates of the corner nodes. For the purpose of illustration, 
the compressible velocity profiles at the Inlet entrance plane, the 
throat and the fan plane are shown in Figures 39 and kO for the QCSEE 
and the Bellmouth inlets respectively. As one observes the radial pro- 
files in Figures 39 and kO the mean flow deviates far from the one 
dimensional model at the entrance plane. The axial velocity component, 
$2 becomes zero at the inlet wall as it is a stagnation point for the 
axial component of the velocity. The radial velocity component, is 
zero at the axis of the inlet because of the axisymmetry of the problem 
and has a large negative value at the inlet wall, indicating the effect 
of the entry lip shape in sucking the flow. At the throat and the fan 
plane is within 5% of QCSEE inlet and within )2% for the 

Bellmouth inlet and also reasonably uniform over the inlet cross 

section indicating that the mean flow is almost one dimensional at the 
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throat and the fan planes. From Tables D.l and D.2 one observes that 
the duct cross sectional area Is minimum at the throat plane for the 
QCSEE inlet and at the fan plane for the Bel I mouth Inlet. Hence, on 
the basis of one dimensional isentropic gas dynamic relations for sub" 
sonic flows one would expect that the axial velocity is maximum at the 
throat plane for the QCSEE inlet and at the fan plane for the Bellmouth 
inlet. One indeed observes in Figures 39 and 40 and Tables D.l and D.2 
that ^2 on the average largest at the throat plane for the QCSEE inlet 
and at the fan plane for the Bellmouth inlet. 

Though the mean velocity profiles at the throat and fan plane 
seem to indicate the presence of almost one dimensional flow at these 
two planes such is not the case at many other axial locations. As 
Tables D.l and D.2 indicate there are very strong radiai gradients in 
and in the vicinity of the Inlet entrance plane and the nose of 
the center body which can not be accounted for by the one dimensional 
isentropic gas dynamic calculations. 

4.6 Acoustic Calculations for the Hard Walled QCSEE Inlet 

With the knowledge of the "compressible" mean flow velocities In 
the Interior of the QCSEE Inlet, one can perform the acoustic calcula- 
tions for the QCSEE inlet configuration. In this section the results 
of the acoustic calculations for the QCSEE Inlet shown in Figure la 
with hard walls will be presented for two frequencies. In order to 
investigate the effects of the two-dimensional nature of the mean flow 

on the acoustic behavior in the inlet, the mean flow was also computed 
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using a one-dimensional, isentropic compressible flow model 


to describe 




Figure 33* Radial Variation of "Compressible" Mean Flow 
Velocities at 3 Axial Locations for the 
QCSEE Inlet. 
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the Inlet meen floM. As In the two dimensional flow computations, a 
free stream Hach nianber of 0.12 and an exit Mach number of 0.52 have 
also been used in the one* dimensional flow computations. The solution 
domain, consisting of the inlet duct, has been subdivided into 3^ 3"node 
triangles with 212 nodes (see Figure 3). The real and imaginary parts 
of the acoustic potential were calculated at the nodes, and the real 
and imaginary parts of the acoustic axial and radial velocity components 
and the acoustic pressure were calculated at the centroids of the tri- 
angular elmnents using the linear interpolation functions for the 
acoustic potential (l.e., see Eq. (52)). 

Calculations for plane velocity wave excitation at the inlet 
exit plane and the frequency u equal 1 and 2 have been performed using 
the l-D and 2-0 steady flow calculations. Figure 41 shows the variation 
of the real and imaginary parts of the acoustic pressure along the upper 
inlet wall, for a frequency of 1.0. The radial variation of the acous- 
tic pressure at the inlet entrance plane, inlet throat and Inlet exit 
plane for u ■ 1.0 are shown in Figure 42, 43, and 44 respectively. Simi- 
lar plots for a plane velocity wave excitation wl th ui ■ 2.0 are shown in 
Figures 45, 46, 47, and 48. 

To isolate the effect of cross sectional area variation alone on 
the wave propagation in the inlet, inlet acoustic calculations for the 
case of no mean flow have been performed for a plane velocity wave excl* 
tat ion at the inlet exit plane for two frequencies, u) ■ I and 2. Fig- 
ures 4g and 50 show the variation of the real and imaginary parts of 
acoustic pressure along the upper wall of the inlet for u * I and 2 
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Figure k\. Acoustic Pressure Distribution Along the 
Inlet Upper Wall (u ■ 1.0; Plane Wave 
Fxci tation) . 


Figure k2. Radial Variation of Acoustic Pressure at the 
Entrance Plane for uj • 1.0 (Symbols defined 
I n Flgur* 1 ) . 


Figure ^3 
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. Radial Variation of Acoustic Pressure at 
the Throat Plane for u « 1.0 (Symbols 
defined in Figure Al). 



Figure kk. Radial Variation of Acoustic Pressure at 
the Inlet Exit Plane for w « i.O (Symbols 
defined in Figure Al). 
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Figure ^*5. Acoustic Pressure Distribution Along the Inlet 
Upper Wall (i** ■ 2.0, Plane Wave Excitation; 
Symbols defined in Figure Al). 




Figure A6. Aadiai Variation of Acoustic Pressure at the 
Entrance Plane for u - 2.0 (Symbols 
defined in Figure A|). 



Figure 47. Radial Variation of Acoustic Pressure at 
the Throat Plane for « = 2.0 (Symbols 
defined in Figure 4l). 
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Figure 48. Radial Variation of Acoustic Pressure at 
the Exit Plane for w ® 2.0 (Symbols 
defined in Figure 41). 



Acoustic Pressure Distribution Along the Upper 
Wall of the Inlet for the Zero Mean Flow Case 
(u) = 1.0; Plane Wave Excitation). 


Acoustic Pressure Distribution Along the Upper 
Wall of the Inlet for the Zero Mean Flow Case 
(for u) = 2.0; Plane Wave Excitation; Symbols 
defined in Figure ^9). 




respectively. Radial variation of the accMistic pressure for the zero 
mean flow case, at the inlet entrance, throat, and exit plane for o) - I 
are shown in Figures 51, 52, and 53 respectively. Similar plots are 
shown in Figures 5^, 55 and 56 for u ■ 2. 

in the problem under coniu derat ion, the propagation of sound 
within the inlet will be affected by the following phenomena: (a) reflec- 

tion from the inlet walls due to Inlet cross sectional area variation; 

(b) reflection at the inlet entrance and exit planes; (c) reflection due 
to the presence of gradients in mean flow properties; (d) convection of 
sound by the mean flow which affects the local effective wavelength; 
and (e) refraction of sound due to the presence of transverse mean flow 
gradients. The presence of some of these effects will be considered 
below by analyzing the predicted inlet acoustic behavior. 

When there is zero mean flow in the inlet, only reflection will 
be important. Furthermore, in this case the Inlet entrance impedance 
* 1 is a good approximation for the almost no reflection condition 
at the entrance plane and thus only reflection from the walls should be 
present. Figures 51"56 show very little distortion in the plane wave 
shape, indicating that there is very little reflection from the inlet 
walls. 

When the mean flow is described by a one dimensional approxima- 
tion, only reflection and convection will be important. Figures 
and i»6-i»8 show very little distortion in the plane wave shape (for 1-D 

mean flow) except near the inlet entrance plane where there are steep 

% 

area and mean axial velocity gradients, reaffirming the conclusion that 


-0.5 


0 


0.5 


1.0 


1.5 


Figure 52 
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. Radial Variation of Acoustic Pressure at the Throat 
Plane for the Zero Mean Flow Case for cu = 1.0 
(Symbols defined in Figure ^9). 




Radial Variation of Acoustic Pressure at the Inlet 
Exit Plane for the Zero Mean Flow Case for u = 1.0 
(Symbols defined in Figure kS) . 


Radial Variation of Acoustic Pressure at the 
Entrance Plane for the Zero Mean Flow Case for 
u) = 2.0 (Symbols defined in Figure k$) . 
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Figure 55« Radial Variation of Acoustic Pressure at the 
Throat Plane for the Zero Mean Flow Case for 
w = 2.0 (Symbols defined in Figure ^9). 



Figure 56. Radial Variation of Acoustic Pressure at the 
Inlet Exit Plane for the Zero Mean Flow Case 
for 0 ) * 2.0 (Symbols defined in Figure 49). 
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there is very little reflection from the Inlet walls. Figures 41, 

49 and 50 show the convective effects of the flow as manifested by the 

overall decrease In effective wavelength for the case with mean flow 

(1-D or 2-D) as compared to the case of no mean flow. Figures 4l and 
45 indicate a minimum effective wavelength near the throat where the 
highest mean flow Mach number exists. 

When the mean flow is described by a two dimensional approximation 
all of the above mentioned effects are present. Furthermore, at the 
inlet entrance plane radial mean flow velocities are quite large and 
there are radial gradients in the axial component of the mean flow. 

Under these conditions the assumed impedance of Z = 1 will result In 
wav.^ reflection at the Inlet entrance plane. For the two dimensional 
mean flow description, Figures 42-44 and 46-48 show severe distortion 
in the initially plane wave shape due to reflection, convection, and 

refraction it Is not possible to isolate the Individual effects. It is 

known, however, that refraction increases with an Increase in fre- 
quency. This effect Is observed by comparing Figures 42-44 (for u * I ) 
and 46-48 (for w - 2) where the distortion of the plane front increases 
with an Increase In frequency. 

4.7 Sound Attenuation Studies for Zero Mean Flow Case 
for the QCSEE, Bellmouth and Cylindrical Inlets 

The objective of this section is to estimate the sensitivity of 
duct attenuation to Inlet curvature and centerbody for a practical tur- 
bofan inlet. The influence of flow gradients on the duct attenuation 
will be included in the next section which deals with the nonzero mean 
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f1(M case. To accomplish this task, a number of acoustic calculations 
on a soft walled QCSEE turbofan inlet, a bellmouth iniet and a straight 
open cylinder of the s»ne overall dimensions (i.e. length and fan diam- 
eter, are performed. The specific acoustic liner impedance values 
chosen are the same for each inlet at a particular frequency and they 
correspond to the near optimum impedance values for a plane pressure 
wave Input into an Infinitely long circular duct. Rice^^ presents these 
values for a wide range of frequencies and liner length to duct diam- 
eter ratios. In the finite element calculations for the zero mean flow 
case a constant acoustic axial velocity is used as the Input condi- 

tion and hence these impedances do not necessarily represent the optimum 
values in the finite element calculations. The radiation impedance con- 
dition is prescribed to be equal to p c. 

Figure 57 shows the dependence of duct: attenuation, ‘^®reduction 

on the frequency of wave propagation, t\{- - -^) for the QCSEE, Bellmouth 

and cylindrical inlets for the case of the full length of the upper wall 

of the inlets being lined. Also shown in the figure are the near optimum 

specific acoustic liner impedances, Z for the four frequencies listed. 

One observes that dB . falls rapidly as frequency increases for 

reouct ion 

all the three inlets which is due to focusing of the sound wave toward 
the inlet axis at higher frequencies. The propagation of high frequency 
plane waves in a duct could be approximated as the propagation in a ray 
tube with its axis coinciding with the inlet axis. As the frequency 
increases the sound power is concentrated near the duct axis and it 
reduces at the wall where it must be absorbed and hence the rapid drop 
In the duct attenuation. One also observes that there is a considerable 
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variation in tiw duct attenuaticMi for the three inlets which can be 
attributed to the differences in their geometries. At lower frequen** 
cies the attenuations differ much more compared to at high frequency and 
fall within 0.5 dB for n ■ 2. in other words attenuation of low fre- 
quency waves is more sensitive to inlet geometry than to high fre- 
quency waves. As explained earlier the scxjnd energy is beamed towards 
the duct axis as frequency increases leaving little sound energy to be 
absorbed at the walls. Hence the propagation of high frequency sound 
waves is insensitive to tl« area variations of the duct resuiting in 
about the same attenuation for all the three inlets. The duct attenua- 
tion for the Bellmouth inlet falls in between the duct attenuation for 
the QCSEE and the cylindrical inlets since the Bei Imouth inlet is geo- 
metrically "median" between the QCSEE inlet and the cylinder. 

k,B Sound Attenuation Studies for Nonzero Mean Flow Case 
for the QCSEE, Bellmouth and Cylindrical InleTs 

The influence of mean flow gradients on duct attenuation is stud- 
ied In this section using the "compressible" mean flow velocities computed 
in Section 4.5 (also see Appendix D). The average fan plane Mach num- 
ber for both the QCSEE and the Bellmouth inlets is 0.32. 

As proposed by Rice,^' the near optimum liner impedance values 

for the cylinder containing a uniform steady flow are calculated by 

2 

dividing the zero mean flow near optlimjm values by (I - M) where M is 
the Mach number of the uniform steady flow in the cylinder. For the pur- 
pose of comparison of the acoustic performance of the Inlets containing 
a mean flow, representative Mach numbers have to be chosen for the 
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cylinder to "represent** in some overall sense the tMO dimensional mean 
flow in the QCSEE and Bellmouth Inlets. For this purpose, one dimen- 
sional isentropic gas dynamic calculations are performed for the QCSEE 
and Bellmouth Inlets and the axial station chosen to *'represent" the 
cylinder is at Z ■ 0.9 for both the QCSEE and Bellmouth Inlets since 
they '*resemble** tiM cylindrical duct at and in the neighbortood of 
this station. The one dimensional Isentropic Mach number at Z ■ 0.9 
for the QCSEE inlet equals 0.579 and for the Bellmouth inlet equals 
0.362. The liner impedance values used in comparing the cylindrical 
duct and the QCSEE inlet are obtained by dividing the impedance values 
listed in Figure 57 by (I - M) where M equals 0.579. Similarly the 
liner impedance values used in comparing the cylindrical duct and the 

Bellmouth inlet are obtained by dividing the impedance values listed 

2 

in Figure 57 by (I - M) where M equals 0.362. in the finite element 
calculations a constant acoustic pressure (p) is used as the input con- 
dition and the radiation condition is again prescribed to be equal to 
p c for the nonzero mean flow case. 

Figure 58 shows the acoustic performance of the QCSEE inlet and 
the cylindrical duct containing mean flow. As in the no mean flow case 
one again observes the rapid fall of attenuation as frequency increases. 
And also the attenuations differ considerably for the two inlets at 
lower frequencies but they lie within 1 dB for n » 2.0. The differences 
in the duct attenuations are due to a combined effect of inlet curvature, 
centerbody and mean flow properties. The duct attenuation for the zero 
mean flow case is higher ccxnpared to the corresponding non zero mean flow 





cas« for n ■ 0.5 and 1.0 for both the QCSEE Inlet and cylinder. For 
n ■ 0. 5 ^'^^reductlon QCSEE Inlet with mean flow Is about 22 dB 

loifi^r than that for the no mean flow case. Also, for n ■ 0.5 
**®reductlon cylinder with mean flow Is about 10 dB lower than 

that for the no mean flow case. However, dB . for the nonzero 

mean flow case for n ■ 2.0 falls within I dB of the zero mean flow case 
for both the QCSEE Inlet and cylinder. These observations indicate 
that though the approximate method proposed by Rice to obtain the near 
optimum Impedance values for the flow case Is quite reasonable for high 
frequency waves, It leads to a rapid drift In the optimum values of the 
Impedances for low frequency waves resulting In a considerable reduc~ 
tlon of duct attenuation. Similar observations are noted in Figure 59 
which shows the acoustic performance of the Bellmouth inlet and the 
cylinder containing mean fl<^. Since the QCSEE and Bellmouth inlets 
are studied for different values of the liner impedance for the mean 
flow case one can not compare their acoustic performance. 

At the present time, the frequency range and transverse mode is 
limited by the maximum number of elements which tne computer can handle. 
In the above cases, seven elements were used to resolve the transverse 
modes and fourteen elements were used to resolve axial variations. For 
a given finite element tr iangul izat Ion scheme, an estimate o'* the accu- 
racy of the numerical solution for various frequencies May be made by 
comparing it with the known analytical solution for the case of a hard 
walled annular cylinder carrying a uniform steady flow. Even though 
there is no one to one correspondence between the accuracies of the 



Bel 1 mouth Inlet 


Cy 1 1 nder 



Figure 59* Oependance of dB^gjj on Frequency for the 
Bellriouth and Cylindrical Inlets for 
Nonzero Mean Flow Case. 
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solution c^tained for a variable area inlet and an annular cylinder the 
above mentioned criterion, in the absense of anything better, is proba- 
bly an acceptable criterion. 


CHAPTER V 


CONCLUSIONS AND RECOMMENDATIONS 
5.1 Conclusions 

An analytical technique utilizing the Finite Element Method 
(FEM) in combination with the Method of Weighted Residuals has been 
developed for predicting the acoustic performance of turbofan inlets 
carrying a subsonic axisymmetric steady flow. An approximate solution 
for the steady inviscid flow field is obtained using an integral method 
for calculating the incompressible potential flow field In the inlet 
with a correction to account for compressibility effects. Acoustic 
properties of the QCSEE inlet, a Bel 1 mouth inlet and a circular cylinder 
for zero mean flow and non-zero mean flow situations have been deter- 
mined for a limited range of frequencies ^2.57) of plane 

wave propagation. Summarizing the results of this investigation it can 
be concluded that: 

1. The finite element technique solutions are in excellent 
agreement with available analytical solutions for the problems of plane 
and spinning wave propagation through a hard walled annular cylinder 
with a constant mean flow. 

2. For an equal number of finite elements, quadrant representa- 
tion of the finite elements is superior to linear representation in 
handling high frequency (to - 15) plane wave propagation. 



3. The duct attenuation and the acoustic pressure distributions 
obtained by the finite element scheme and the finite difference results 
of Baumeister for the case of a lined rectangular duct carrying a 
steady uniform flow are found to be generally in good agreement. 

k. The results obtained by the finite element program using the 
"exact" Impedance at the open end of the duct provided by the integral 
solution approach of Bell, Meyer and Zinn are found to be in excellent 
agreement with the results of the integral solution approach for plane 
wave propagation in hard and soft walled QCSEE inlet carrying no mean 
flow for a range of frequencies. The simple "no ref lect ionf' impedance 
condition at the open end of the inlet though inaccurate for low fre~ 
quencles {id ^ 2) of plane wave propagation, approaches the "exact" impe 
dance condition for high frequencies (u) = lO). 

5. Results for low frequency (u ^ 2) plane wave propagation 
through the hard walled QCSEE inlet containing a mean flow show that 
when one- dimensional steady flow is assumed to exist In the inlet, the 
plane wave propagates with relatively little distortion. However, pro- 
pagation of a plane wave through the fully two-dimensional flow field 
in the inlet produces severe distortions due to the excitation of 
higher order modes. 

6. Plane wave calculations for soft walled QCSEE inlet, a cir- 
cular cylinder and a Bellmouth inlet for the near optimum liner impe- 
dance values of an infinitely long circular duct with no mean flow, 
indicate that the duct attenuation falls rapidly with increase in 
frequency for the cases of zero mean flow and fully two-dimensional 
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axisynvnetric mean flow due to focusing of the sound wave toward the 
duct axis for high frequencies. Attenuation of low frequency plane 
waves is found to be tnore sensitive to inlet curvature, center body and 
mean flow gradients as compared to that of high frequency plane waves. 

The approximate method proposed by Rice to obtain the near optimum liner 
impedance values for the flow case from the zero mean flow case Is 
found to be reasonable for high frequency plane waves but leads to a 
rapid drift in the optimum values for the low frequency plane waves 
resulting in a ccxisiderable reduction of duct attenuation. 

5.2 Recomittendat ions for Future Research 
Several improvements over the developed theoretical analysis 
which could be made to broaden its range of applicability are discussed 
below. 

I. Prediction of the Far Field Noise Levels 

The present analysis predicts the ratio of the input acoustic 
power at the fan plane to the output acoustic power at the inlet entrance 
plane and thereby concludes about the effectiveness of the liner to 
absorb sound. However, in practice, the effectiveness of the liner is 
estimated by measuring the sound pressure levels in the far field. 

The far field predictions could be made using the finite element 
technique by extending the tr iangu 1 izat ion scheme from inside the inlet 
to the "far field" through the inlet entrance plane. The elements 
external to the inlet could be made larger as one moves away from the 
inlet since the variation of the acoustic variables in the open space 
is not expected to be as large as it is inside the inlet. In this 



scheme the rad tat ton condition should be prescribed at the peripheral 
elements In the open space. The prescription of the simple p c radia- 
tion condition out in the "far field" Is expected to be quite good 
since any wave front approaches a plane wavefront locally as one moves 
sufficiently far away from the source. The extent of the finite ele- 
ment spread out in the open space to simulate the far field will be 
governed by the ccxnputat ional capabilities and mass storage space of the 
available computer system. 

il. Effect of the Molecular Transport Properties on the Acoustic 
Behavior of the Inlet 

In the present theoretical model the effect of viscosity and ther- 
mal conductivity of the fluid on both the mean flow field and the acous- 
tic field have been neglected. These assumptions and the assumption of 
irrotational flow simplified the boundary value problem to the solution 
of a single complex partial differential equation with variable coeffi- 
cients for the acoustic velocity potential subject to a set of complex 
mixed type boundary conditions. When employing the finite element tech- 
nique the use of the velocity potential offers many advantages over the 
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conventional linearized gas equations approach. For two dimensional 
mean flows the velocity potential approach reduces the computer storage 
and running times by an order of magnitude compared to the more general 
linearized jas equations approach. 

The acoustic equations for sheared viscous flows were developed 
33 

by Mungur et al. Since the gradients in the wave motion are "small" 
the molecular transport of momentum and energy due to viscosity and 
thermal conductivity of the fluid respectively are neglected by common 
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practice in deriving the linearized acoustic equations, and the inviscid 
acoustic equations are solved using the mean flow velocity field that 
results from a consideration of viscosity in the mean flow field only. 

Because of the significant computational advantage of the acoustic 
potential formulation, one is strongly motivated to extend the acoustic 
potential finite element program to include the case of sheared mean 
flows. Goldstein and Rice suggest a very practical procedure for cal- 
culating the effect of the boundary layer in conjunction with the poten- 
tial function analysis. Goldstein and Rice point out that in many cases 
a linear shear layer can be used to accurately model a 1/7^^ power law 
turbulent boundary layer. For the special case of a linear shear bound" 
ary layer, they show that the following correction procedure can be 


app lied. 

The specific acoustic liner impedance, at the inlet wall is a 
known parameter. For a given boundary layer thickness 6 and a mean flow 
Mach number at the edge of the boundary layer, a specific acoustic impe- 


dance at the edge of the boundary layer, 


(Z„) rr can be calculated 

£ effective 


which is independent of'the potential core flow. This effective impe- 


dance could be used as the boundary condition for the potential flow 


finite element calculations. The procedure is approximate because the 

34 

correction procedure assumes a uniform boundary layer in the absence 


of mean flow variations. Nevertheless the procedure represents a simple 
and effective way of accounting for shear with a potential approach for 
the duct acoustics problem. 
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III. Optimization Studies for the Liner Design 

The present technique predicts the **®reductlon ® *^*J®*^ 

containing a two dimensional axisymmetric mean flow for a prescribed 

conf iguraticMi of the sound absorbing liner. The liner performance 

which could be measured as the dB_ . . . per unit length of the liner 

reduct i on 

is a function of its specific acoustic impedance, frequency of wave 
propagation, mean fltx^ characteristics and the geometry of the duct. 

The technology needed for a practical liner design has to predict the 
optimum liner configuration for a desired Such a design 

information could be prepared by conducting optimizaticxi studies for 
liners over a realistic range of the independent variables. Such a 
study will provide the needed technology data base for the liner 
design. 
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APPENDIX A 


PROOF OF EQUATION (60) 


To prove 


1-// at)^9'>(L*)VdZ.-j^-y^x2. 
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Consider the linear transformation 
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The elemental area drdz is related to dL.dLj through 
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where the Jacobian J 
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Substituting Eqs. (51) in (A-3) one obtains 


L(e) ^_{e) 

1 j 


(z, -zkHr.-r,^)- (r,-r,^)(z.-z^) 


(A-k) 


2x (Area of triangle ijk) * 2A 


118 


Hence, Eq. (6o) becomes 
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That is, 
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Using the definition of Beta function and Ganma function, the inner 
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integral is equal to 
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Again using the same formula above, 
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APPENDIX B 

DEVELOPMENT OF ELEMENTAL REUTIONS 

B. 1 Application of Galerkin's Technique to a Finite Element 
Equation (59) when applied to a single finite element yields 

III N^^®^L 4 ')dV^®^ - 0 , p - I, 2 ,...,o (B- 1 ) 



where a is the total number of nodes in the element e (o is 3 for a 

linear and 6 for a quadratic element). For illustration purposes the 

Galerkin's technique Is applied to Eq. ( 6 ) in this development. Exactly 

same procedure applies to Eq. (7) also. Hereafter, the superscript e 

is dropped since this derivation is for one el&nent only. 

The first three terms in Eq. ( 6 ) contain second order partial 

derivatives of i|) which are reduced to first order by using Green's 

theorem in a plane. Green's theorem in a plane region R bounded by a 

24 

curve C for two functions U(r,Z) and V(r,Z) is (see Figure B. 1 and note 
t^at C is traversed in a direction that R lies to the left of C) 

U [I 7 - -^] drdz » I Udz + S Vdr (B-2) 

R 3z 3r 

Ident i fy ing 


and 
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V - N^r^[(c^- 

and applying Eq. (B-2) yields after some manipulations 

- - - - ll - - 2 

^ ^^rzh ’ T" * (Y-D 7 -+ (Y-l)#^« 2 ^}N^r^] 

♦ *2% ■ ♦!» 

+ ((v+l)i^^ij ♦ * (Y-OJrr*z ■ 

2-2 

- - ^2-5- ]N + 2wl N r^ i + 2wi N r^ i 

P 'M ' 2 U 2 

+ u(Y-l)[i^j. + 7- + 

- f N^r^Uc^ - ^1)^2 - 2i^i^^^]dr-^ " if)0^dz, 

C 0 

U“I,2,...,cJ ( B“ 3 ) 

It is simpler to evaluate the boundary integrals in Eq. (B-3) in terms 
of the natural coordinates (s, n) where s Is the coordinate along the 
curve C and n is the coordinate normal to C and pointing away from the 
region R (see Figure B-’2). Differential lengths dr and dz along the 
curve C are related to ds by 

dr ■ ds sin a 

(B-M 


dz « ds cos a 
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and the velocity conn>onents are related by 

0 * 0 sina ~ 0 cosa 

s " (B-5) 

0 ^ “ sJ®® 

Substituting Eqs. and (B~S) in the boundary integral of £q. (B-3) 

one obtains 


^ ” ^2)^050 sina - 20^0^ sln^al 

M ^ W S 


0 ^ {c^ - (0^sina - 0^ cosa) }]ds , w “ 1, 2,..., a (B“6) 


it is to be noted that the contributions to the boundary integral 
in Eq. (B-6) come only from the elements located on the boundary 
surfaces of the inlet as the individual element contributions located in 
the interior of the inlet are cancelled as the element boundaries are 
traversed in exactly opposite directions for the two neighboring elements. 
Hence, though the individual element boundary integrals will be nonzero, 
when the elements are assembled into a global matrix form, the internal 
element contributiotiS to the boundary integrals vanish. So Eq. (B-6) 

Is evaluated only for elements located along the boundary of the inlet. 

B.2 Element B e I ations for an Internal Element 
For an internal element the right hand side of Eq. (B-3) Is not 
evaluated as discussed above. The left hand side of Eq. (B-3) is eval- 
uated by using linear interpolat Ion functions for 3~node triangles and 
quadratic interpolation functions for 6-node triangles (see Eqs. 51"58). 
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B.2.1 Evaluation for a 3*Node Triangle (See Figure 5) 

For a 3“node triangle the variation of the acoustic potential 
over the element is given by 

^(r,z) - L.(r,z)|j + U(r,z)^^ + Lj^(r,z)|j^ (53) 

and 

|(r,z) - L.(r,z)$. + Lj(r,z)|j + Lj^(r,z)||^ (5^) 

For the sake of simplicity in evaluating terms in Eq. (6-3) the mean 
flow quantities c, and v/ithin each element are assumed to be 
constant and equal to the aritNnetic average of their respective values 
at the nodes i, j and k. However, the derivatives of the mean flow 
quantities like c^, c^ etc. are evaluated assuming that 

the mean flow quantities are varying linearl/ over the element just like 
the acoustic potential. Also the coordinate r is assumed to have a 
constant value r equal to the r-coordinate of the centroid of the element 
except where the derivatives with respect to r are to be evaluated. 

The following polyncxnial integrations of the area coordinates need 
to be performed: 


a. // L drdz, p 

A 


i , J and k 


and 


b. jj + Lj(}ij + Lj^<f,j^]drdz, p = i , J and k 


Application of Eq. (60) to integrals a and b yields 


12^ 


// L drdz - ^ , 


and 


u = 1 , j and k 
n = i , j and k 


0 A 


y ri j'j 


1 '^j • '^k ■ ''un'^n 


(B-7) 


where 6^^ Is the Kroneker's delta having the property 


yn 


•! 


1 if y = n 

0 if y ri 


(B-8) 


Substituting Eqs. (B-7) in Eq. (B-3) one obtains the desired element 
relation for an internal element: 


(c.$. +c.i. +Cj^^l^) fc 

i^4 X 


Pu -2 ,-2 22. -2 -2. 

^ r (c -<j.^)+^r(c -4,^) 

_ ^ w2 

+ r (c Cj, - ~ ^ 

r 

- - r^ X 2A ~2 - - 

+ (y- 1)^+ i-p^-2rM-^<^^^^ 

r . . 2 


[- 


^ ^ 4, ) ¥> + — ^^“1^ 

^rz^z ^r^zz 6 2A 


(c - 4>2<I>2 z) ( (y+ 1 )4>224’z'‘‘2^rz‘^r 


+ (Y-l)i^^^2 ■ 


zV, 2A 

— ) ’'T 

r J 


- ^ )‘r^ X li 

+ 2o)4,^r^ !Ui- + Zco^^r — ^ 

+ u)(Y-l) X [^^^ + -L+ (^. +;. + X 0, 

y = i,j and k n * i.j and k (B- 9 ) 
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B.2.2 Evaluation for a 6- Node Triangle (See Figure 5) 

For a 6-node triangle the variation of acoustic potential over 
the element is given by 


(|)(r,z) - N. (r,z)i. + Nj (r,z)(|)j + N|^(r,z)i|^ + Nj^(r,z)|j 


+ N (r,z)4> + N (r,z)(|> 

m m n n 


(55) 


and 


$(r,z) = N.(r,z)^j + Nj(r,z)^j Nj^(r,z)^l^ + Nj^(r,z)^^ 


+ N (r,z)if> + N (r,z)<() 
m m n n 


(56) 


where for corner nodes 


N, = 2L^ - L, 


N. 

J 


= 2LT - 
J 

2 


L. 

J 




(57) 


and for mid- size nodes 


N„ = i»L.L. 

^ I J 


N = 1*L.L, 
m J k 


(58) 


N = kl, L. 
n k I 


The same assumptions regarding the mean flow quantities and the 
coordinate r made in Section B.2.1 are applied to the quadratic elements 
also. 

The following polynomial integrations of the area coordinates 


need to be performed: 
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I “ // N 4 drdz 
ay V y r 

a 

I. “ // N (j) drdz 
by V y^z 

A 

3N 

I “ // 7 — ^ 4* drdz 

A 

u = l,j,k,£,m and n (B-IO) 

'du ■ // ^ *r 

A 

3N 

I = // T— ^ 4> droz 

ey ■'/ 3z ^z 

A 

•fy = // % ^ tlrdz 


The integrations will be performed for y = i and y = £.. The integra- 
tions for y = j and k and y = m and n can be obtained by a proper cyclic 
rotation of indices i, j and k and £, m and n (see Eqs. (52))- 
Case 1 y = i 

Cons ider 


I = // N. (J) drdz 
a. •'/ I ^r 
I A 


Substituting Eqs. (55/56), (57) and (58) into I yields 

a . 

I 


a . 

I 


I, = //(2Lt - L,) X H 


(l|L,-l)c.d.+ (<lL.-nCjtj+ 


+ ^(c.L. + c,L.)(J). + A(c.L, +c,L.)<j) 
ij Ji S, jk kj^m 


+ 4 (c. L. + c. L, )(}) I drdz 
k I I k 


n 
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simplifying and evaluating the Integrals using Eq. (60) yields 

1 , “ "57: t2Cj4», - c. (j); - c, A, + A (2c. - c. ) - A (c. + c, ) 

at 30 rl j^j k^k j 1 ^m j k 

+ 4>n(2c^- c.)] (B-11) 

Consider 


1. . = // N. A drdz 
bl V I 
A 


Substituting Eqs. (55/56)» (57) and ( 58 ) into I. yields 

b • 

I 

+ (bL,^- l)b|^*l^ + b(b.L. + b.L.)4,^ b(b.L|^+b|^L.)»^ 

+ (b, L. + b. L, )(|) ] drdz 

k I 1 k ^n 


Simplifying and evaluating the integrals using Eq. ( 60 ) yields 


I. . - "TTr [2b. 4>. ~ b.(J). ~ b,(j), + (2b. ~ b.)i}) - (b. + b, )a 

bl 30 1^1 yj k^k j i j k 




(B-12) 


Cons ider 


I . 
ci 


9N. 

/ h — ^ <t> drdz 
■'/or ^ r 

a 


Substituting as above from Eqs. (55/56), (57) and (58) into 1^. yields 
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* "'"■I" * ('"■j-')Cj«j 

+ %- DVk 4 (c,u + c.l,)»^ + 

* ‘‘<‘=kS * 'i'-k>*n' 


SimpUfylng and evaluating the integrals using Eq. (60) yields 



Cons ider 


3N, 


di 


II -5r V 


Substituting as above from Eqs. (55/56), (57) and (58) into I .. 

d I 




+ (^L, - l)c,d», + 4(c.L.+c.L.)(b + ^(c.L +c,l 

k k^k ij Ji £ jk k 


+ 4(c.L. + c.L, )di ] drdz 

k I I k ^n 


Simplifying and evaluating the integrals using Eq. (60) yields 

'di ■ Ia * * Vk'" i’ 





(B-13) 


yields 



(B-lM 


129 


Consider 

3 N 

I .= ff 6 drdz 
ei V 8z '^2 
A 

Substituting Eqs. (55/56). (57) and (58) into 1^. yields 

I > //— W (H - Dbix [(4L - 1)b * + (4 l - l)b 
ei ^ ^2A)^ ' » I • J J J 

+ ('"■k- l)b|^*|^ + 4(b,L. + b.L.)*^ + 4(bjL|^l-b|^L.)4,^ 

+ 4 (b, L, + b.L, )(j) 1 drdz 
k i I k ^n 


Simplifying and evaluating the integrals using Eq. (60) yields 


I . 
ei 


= X lb.*, (i) + bj*j (- + b|^*|^(- i) 


+ b.(f) 

ysi 


(f) 


+ 4 > (O) + b 


' m 


k^n 



(B-15) 


Consider 


I.:. = // N. (j) drdz 

A 

Substituting Eqs. (55/56), (57) and (58) into t|:. yields 

Ir. = f/(2L? - L.) X [(2L? - L.)(j). + (2L? - L.)<(,. 

f I A ' I I I J J J 

+ (2lf - L, )(j), + 4L.L.d,^ + 4 L.L. 4 . + 4L, L.a ] drdz 

' k k^k t J k*m k i^n 

Simplifying and evaluating the integrals using Eq. (60) yields 
If. = 24 X + ^j(- 3 I 0 ) + 3 ^) + ^j,(0) 


(B-16) 
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Substituting Eqs. (B-li) - (B-l 6 ) into Eq. (B~3) yields the following 
quadratic element equation for an internal element for the corner node i 

- . 2 - 2 -2 - - - - *2" - _o- - 

I , X [2r (c* ” 4 ) + 2r‘ (cc “44 )“2r4 4 , - 2r4 4 ,, 

ai ^r r ^r^rr ^rz^z zz 

-2 4^ 2 

+{(Y+l)4rr*r ^^rz^z ' ^ * (Y“ 1 )^+ (y - 1 )4^422 >»■ 1 

r r 

*'bi * K ' ♦z^zz* * 

-2 

■j ~2r 2 2 c 1 . i o “2 T i '1 o ~2 ~ 

- I^.r [o) - m ry -I Igj 2a) r 4 ^ + 1^^. 2o) r 4 ^ 

r 

9 ^ 

+ 1^. oj(Y-l)r [ 4 ^^ + -p + 4 ^^] = 0 (B- 17 ) 

where the superscripts - and ■' refer to the real and imaginary parts of 
the complex integrals in Eqs. (B-11) - (B-I 6 ) respectively. Element 
equations similar to Eq. (B-17) for nodes j and k may be obtained by the 
cyclic rotation of indices i, j and k and I, m and n in the Eqs. (B-11) 
(B-16). 

Case 2 p = & 

Following exactly the procedure indicated in Case 1 one obtains 
the following expressions for the integrals in Eq. (B-lO) for y = l: 
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II tlrdz 


+ ^ {Cj + 2C|,)4,^ + tI (<=i * 2‘k>*n’ 


(B-18) 


•bt * // % *2 ‘"■'‘^ 

A 


K . ^ 


[b|»,(To> + »j*j W * Vk<- 30> * Ts * ^•’♦b 

" TT <‘'J " 2'>k>*m " iT “’! * 


(B-19) 


I « // -r— ^ 4> drdz 

cH V 3r ^r 
A 

- 2? '*i * *j " ♦k*®) ^ T <‘^i * ^i-^j " 

+ X ■* ^'"i'=k *'=■]* * T 


+ Zc.c^ * c.c,)l 


9N 

'dt = (/ ^ 

-2f’‘ l*i(^)"*j<^>**k<®> *‘■i"j"^■'=i 


+ 2bjC.) + -f (b.c. + 2b, C|^ + b.c. + b.c,^) 

X <'>i\ " 


(B-21) 
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3N, 


I . “ // 6 drdz 

ei ■'/ 92 ^2 


+ (b.b. + 2b. b. + b. + b.b,) + (b . b. + b , 

o lj ik j jk 6 ik 1 


+ 2b. b. + b.b.)] 

J k J I 


(B-22) 


•fji " //Np (j>drd2 


2A [<J).(0) + (|)j (0) + -^) + 




(B-23) 


Substftuting Eqs. (B-J8) - (B-23) into Eq. (B-3) yields the following 
quadratic element relation for an internal element for the mid-side node 

X [2?(c^ - ^J) + 2r^(c - ^r^rr^ “ ^rz^z ' ^’’r ^r^zz 

-2 

- . '2 d,^ . . ,2 

+ { (y+ 1 )^rr^r * ^^rz'^z " ^ J 

r r 

+'*b£ t2r^(c ^^rzh 

+ (Y-l)irr^z - (Y-1)^} r^) 

r 

■^^c£ ’ ^d£ V^z^ *e£ 
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-T,, ^ 2„-r^ 

+?^j^ w(r-l)r " ® (B-24) 


where the superscripts ” and - refer to the real and Imaginary parts 
of the complex Integrals In Eqs. (B-18) - (B-23) respectively. Element 
relations similar to Eq. (B~24) for mid~side nodes m and n may be 
obtained by the cyclic rotation of indices i, j and k and m and n 
in the Eqs. (B-l8) - (B-23). 

B.3 Boundary Integral Evaluation for a Hard Wall 
The boundary condition at a hard wall of the duct in terms of 
real and imaginary parts is 

(fi * 0 

and (15) 

A » 0 
n 

B. 3.1 Evaluation of Hard Wall Boundary Condition for a 3 "Node Triangle 
(See Figure B.3) 

Substituting Eq. (15) in the boundary integral expression given 
by Eq. (B-6) one obtains 


8 





- i^) 


cosa sina " 2(|i d» sin 


r’^z 


’a] ds 


(B- 25 ) 


For the 3"node element configuration shown in Figure B-3 the interpola- 
tion functions along c reduce to 


13 ^ 


Nj - 1 - 

Nj - C [0 1 ^ 1 1 (B‘26) 

N, iO 1 

Hence the variation of <;> along c becomes 

^ - (I “ Oi. C0j (B-27) 

Application of the chain rule for differentiation yields 

igds - {^j - (J.j)d^ (B-28) 

Substituting Eq. (B-28) into Eq. (B-25) one obtains 

"■ — 7 ^ 

6j ■ r [{<!> - ((> )cosa sina - 24 > 6 ^ sin a] / (l-C)(i.-i.]d^ 

j " ^ ^ 5-0 (5) J ' 

That is 

- “2 -2 “2 - - 2 i ~ ^ i ^ 

S - r [(4> - (}) )cosa sina - 2<|! 0 sin a] — ■ — , 

M r z ^r z 2 

p = i and j (B-29) 

Since and are unknown as of now, 3^ is transported to the left 
hand side of Eqs. (B-3). 

B. 3.2 Evaluation of Hard Wall boundary Condition for a 6 -Node Triangle 
(See Figure B-3) 

The boundary integral to be evaluated Is given by Eq. (B-25). 

For the 6-node element configuration shown in Figure B-3 the interpola- 
tion functions along c reduce to 
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Hj - 25 - 3C + 1 


N. - K 


M i. 5 1 * 


N. - N - N =0 -* 

k m n 


(B-30) 


Hence the variation of $ along c becomes 

i - (2c^ - 3e + D^j + (2C^ - C)<!>j + ^(C - (B-31) 


Application of the chain rule for differentiation yields 


- [(^ - 3)$j + - 1)^5 ♦ 1*(1 - 2c)^JdS (B-32) 

Substituting Eq, (B“32) into Eq. {B-25) and evaluating the polynomial 
integrals leads to 


”2 *2 ”2 “ “ 2 
3 , ■ r* [ (4»^ * ‘I’z ^ cosa sina " ^ 

lij (" j) * (’ ij (j)3 


_ ” 2 “ 2”2 .2 
3 j = r [( 41 ^ " (ji^)cosa sina " 4>2 ^ 




(B-33) 


(B-34) 


and 


-2 


•"£*2 ^ 
r^ [( 4 ^ " 4>2 ^ cosa sina " 21^^42 * 


II j(” y) + Ij (-j) + 


(B-35) 
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As i|, and 0 ^ are un nown, B., Sj and 3 ^ are transported to the 
left hand sides of the corresponding equations in (B~3) for nodes i, j 
and i respectively. 


B.k Boundary IntfcMral Evaluation for a Soft Wall 
(See Figure B-3) 

The boundary conditions at a soft wait using the concept of con- 
tinuity of particle displacement are given by 




. . - . P *5 - - - 

“p[mi + <fc i ] + " 0 + A 4 1 

s^s U) s ss s s ss 


* t It IT 


(29) 


and 


p4i ^ __ 

X(fe +0* " “P (1) i 1 “ — - [u)$ + ^ 1^) + 41 ^ 1 

n ^n s’^s u) s ss s s ss 


. (^2 ) P. i£ {u,; + $ ^ ] 
>-r - 9s w ‘ ^ 


(30) 


In the absence of mean flow the above equations reduce to 



(B-36) 

~ 

(B-37) 


Since Eqs. (B- 36 ) and (B- 37 ) do not contain the second derivatives of 
the acoustic potential, one could use 3"node triangles to evaluate the 
soft wall boundary conditions. 
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B.4. 1 Evaluation of Soft Wall Boundary Condition for a 3‘’Node Triangle 
(Nc mean flow case only) 

The boundary integral given by Eq. (B-6) reduce to the following 
for the case of no mean flow: 


B 


y 



(B-38) 


The soft wall boundary condition for the case of no mean flow given by 
Eqs. (B- 36 ) and (B-37) are solved for (ji^ and and on substitution into 
Eq. (B- 38 ) leads to 


3 = 

y 


I N r 
C ^ 


2 


^ 2 
Q + X 


Ids 


(B-39) 


Using the interpolation functions given by Eq. (B"26) for the element 
configuration shown in Figure B-3 one may write 

^ = (1 - 

and (B-1»0) 

$ = (1 - 

The differential length ds along C is given by 


ds = ((Zj - Z.)^ + (r^ - r.)^l’^^d^ = L.jde (B-4l) 


Substituting Eqs. (B-AO) and (B-^1) into Eq. (B-39) yields 


e = / 

y ^ 


-2.. 

r N 


^-0 


•U)0{ (l-^)^j + (l-c)<^. + ) 

2 2 
0 + X 


L.jdE 


M = i and j 


Substituting for N from Eq. (B-26) and evaluating the polynomial 


integrals yields 


6 {0^ + x^) ^ 


-e($. + «^. + 6 I ) 


+x(^. + A. + 6 ^ ) 

yn n 


p = i and j n » I and j 


Since the nodal values of the acoustic potentials are unknown, 3^ is 
transported to the left hand side of Eq. (B-3)* 

B.4.2 Evaluation of Soft Wall Boundary Condition for a 6-Node Triangle 
Since the soft wall boundary conditions given by Eqs. (29) and 
( 30 ) contain second derivatives of acoustic potential in the presence 
of a slip flow at the liner, quadratic elements or higher order ele- 
ments are to be employed to treat such boundary conditions, Eq. (29) 

and (30) can be solved simul iianeously for A and A to yield 

^n n ' 


n ^2 /. 

C 1 .< 


y [v- X {- p(u$ + + 


(-wi + i ^ + J ^ ) 


ss's s ss 


It TT + ^ $ ) } 

Y I ^ dS U) S S 


+ X X {- p (“u)A + (j) A ) 
S S 


(w^ + (j) ^ + (j) ^ ) 

s ss s s ss 




(B-43) 


I 
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- ii ~T ^ 


-X X {-p(4 + (.^1^ + iss^s'"is^ss^ 


Substituting Eqs. (B~ 32 ), (B~l*l) and (B~l< 3 ) into the boundary inte~ 
gral in Eq. (B-6) and regrouping the terms one obtains 


-w — 2 *“ 2 "■ 2 

* *" ^ (‘i’r ~ <{>7)^050 sina - 24, a- sin^oil / N * dr 

C=0 ^ ^ 

-2 

r “2 “ - ?] 

+ 2 2 ~ ~ 4 ’f ^ij ^ 




*s*ss , 2 , ^ ac X -2 
♦s* 2 P *56 ;;; x - (yrr) - iF “ ♦ > 

c ^ s 


A - - p ^c^cc 7 ~ '^r n -7 

+ 4>g{“2p 4>e^ ^ 6 + ^T~T^ fi- i- i 4>^J 

5 s 0, '“I : 8s to S-' 


(^ly) ^ If 4 > s 9 


+ 4i{“(y_|) ^ gg <l>5^ " pa)6}] ^ = 1 »j and £ (B-l»5) 



The required polynomial integrals in Eq. (B-45) are evaluated below 
by using Eqs. (B“30). 


I 1 , , 

“ / N.d? = / {2C - 3C + l)dC » 7 

C-0 ' E=0 ^ 


I I 5 , 

S. = / NdC-/ 2C^-CdC = l 
J C=0 ^ 


I I - 

^a£ “ ! " / ^Ul-5)dC - 

^ 5=0 ^ E=0 ^ 


S. . = / N.(}. d = / (2C - 35+ 1)[(4C- 3)0. + 1)<^ 

' 5=0 ' ^ ^ 5“0 ' 

+ 4(1 - 25)0^]d5 
= 0.(- j) + 0j(- |) + 0^(f) 


S . =/ N.0,d5 = / (25^ - 5)I(^5- 3)0. + (^5 - 1)0. 

bj ^=0 J ^ 5=0 ' J 

+ 4(1 - 25)0j^]d5 

= 0; ij) + 0j(j) + 0/- |) 


1 1 . 

Sko = / N.0^d5 = / 4(5-5 )[(45-3)0. + (^5- 1)0; 

5=0 ^5=0 ■' 

+ 4(1 - 25)0^]d5 

= n (-f)+0j(f) + 0,(o) 

S =/ N.0d5 = / (25^- 3? + 1)[ (25^ -35+ 1)0;+ (25^- 

5=0 ' 5=0 


(B-46) 

(B-47) 

(B-48) 

I 

(B-49) 

(B-50) 

(B-51) 

5)0j 


\ti\ 


s . 

cj 


1 

/ N <{.de 

^-0 J 


/ {2C^ - OI(2C^- 3C+ + (2C^- C)4>: 

4-0 ' 


+ 4(4 - 4 )(^^]d4 


<l>j {- 3 q) + ( 15 ) + 


(B-53) 


C£ 


/ Nji'}>dC 

C=o 


/ 1,(4 - 4^)((2C^- 3^+0*. + 

4=0 ' J 

+ 4(C-4^)4.Jd? 


,(ii> * ♦j(-ri) - *.( 1 !) 


(B-54) 


Also from Eq. (B- 30 ) one may obtain 


) = -= — ['}>. + ‘l>. 

ss |_2 I j 

ij 


2^p] 


(B-55) 


Substitution of Eqs. (B-46), (B-49) and (B-52) into Eq. (B-45). for 
u = i yields the desired evaluated boundary integral 


"*2 *2 *2 - - 2 * 
Si = r - 4^)cosa sina - sin a] S 


bl 


+ • -2 - - 2 
— [c^ - (d) sincx - 4> cosa) ] x 

2 2 r 


0 +x 


[*ss<- 


- X . + i (-0 i^)s . 

u) s ai ss CO s a I 


+{-2p (1)^0 


s ss 


, 2 V p 8c X -2 , bi 
X - <FT> r 3i- - *5 > — 

C IJ 


p 4> 

. r oT 7 . S SS 

+ I -2p (}i X + 

s to 


e .. (2 

T - 1 - 8s (0 s . 

C IJ 


\k2 


i,e^;u.x) s^, 

0 

+ (‘( 77 ]*) ^ci ^] 


where the superscripts - and '• refer to the real and imaginary parts 
of Eqs. (B-49) and (B-52) respectively. 

Similarly, substitution of Eqs. (B-l»7) {B“50) and (B-53) into Eq. 
(B-^5) for yields 


- ~2 -2 “2 “ “ 2 ~ 

3j * f sina - sin a] Sj^j 

-2 - - 2 

+ 2 ~ (4*2 “ ^r ^ 

e + X 


L^ss 03 aj ^ss 0 ) aj 


^ '^s'^ss 

+ {-2p *0 - ^ 

s 03 


- - P '*’s‘^ss 

+ {-2p ^ X + ^ 

3 (i) 


^ W-r - 3S 03 U 

C I J 

a + ( ■--) ^ 1 

® W-r - 3S 03 '*“s ^ L,, 


^ 3S 03 - -Jj 


+ { - (yrj’ - if *5® ■" 

C 


+ {-(--f) f is"" " 


(B-57) 


On similar lines, substitution of Eqs. (B-^8), (B-51) and (B~5^) into 
Eq. (B-45) for p = £ yields 
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((ij - Cosa sina - 

r\.. 

+ 2 * ■^2 " ^r » 

e + X 

[a (- £• X i^)s + i (^e i^)s 

^^ss ' U) a£ ^ss^ui ^ s' al 


. f T ^ ^ *^5^SS / 2 X p 3c X -:2^ ^b£ 

+ {- 2 p ^ e — - — X - (tt> - ir ^ * 5 ^ rr 

C IJ 

. f «" 7 . ^ ^s^ss „ . / 2 ^ p 3c 0 72 t ^b£ 

+ i-2p (f) X + 0 + { — rJ ~ T / 1 — 

^ w 'y- 1 - 9s CO L, . 

c !j 


^-{-(^rT>r|f + p“x) 

c 


+ ^ ^ (B-58) 

C 

As 4>., <J>. and <|). are unknowns, 3., 6. and 3„ are transported to the 

I J 36 3 J ^ 

left hand sides of the cc rre spend ing equations in (B-3) for nodes i, j 
and ^ respectively. 


B. 5 Boundary Integral Evaluatior *^o ‘ Fan Plane Source Condition 
The sound excitation condition at the inlet fan plane could be pre- 
scribed as velocity source (Eqs. 10 and 11) or a pressure source (fqs. 

13 and 14). The boundary integral will be evaluated for both types of 


source conditions. 


\kk 


B.5»1 Evaluation of a Velocity Source Fan Plane Condition for a 3~Node 
Triangle (See Figure B.4) 

Along the curve C a Is equal to 90** and assuming that the radial 
component of the mean velocity is zero (which has been found to be 
reasonable from the mean flow computations) the boundary integra'. In 
Eq. (B-6) reduces to 


6 

P 





(B-59) 


The boundary condition in terms of a radially varying velocity source 
is 


and 



|,»f(r) 


( 10 ) 

( 11 ) 


since the direction of the outward normal is in the positive z direc 
tion. Using the interpolation functions given by Eq. (B-26) one may 
wr i te 


♦n - 

r = (l-C)r. + Cr^ ^ ^ ^ ^ ' (B~60) 


ds = (rj - r.)d£. 


Substituting Eqs. (B“60) into Eq. (B-59) one obtains 




Evaluation of the various polynomials in the above integrals yield 
the following forcing boundary conditions for the linear element in 
terms of a velocity source which form the right hand side of Eq. (B-3) 


and 


8| - X [f,{r|(i+ r,_r. (^)trj(^)) 


* V^i*20* * '' 10 * 15 * * ''j*20**' 




+ fj('-|(35) + r.r.(fi)+r?(i))] 


B.5'2 Evaluation of a Velocity Source Fan Plane Condition for a 
6-Node Triangle (See Figure B-4) 

The boundary integral to be evaluated is Eq. (B-59)- Making 

of the interpolation functions for quadratic elements given by Eq 

(B- 30 ) one may write 

= (2C^ - 3c + Dfj + -Ofj + ^(c- 

r = (1 -e)r. * ir. 0 ^ 

ds = (r . - r . )d^ 

J ‘ 


(B-61) 

(B-62) 

use 


1 1 
(B-63) 


Substituting Eqs. (B-63) into Eq. (B-59) one obtains 



1^6 


6| 

j 

I 


1 



(2c^- 3c+ 1) 
4 ( 5 . 52 ) 


[ {(l“5)rj + ir. f(c^ - X 
(25^- 35+ l)fj + (2C^-C)f^ + i*(C- 5^)f^)l 


X (pj ‘ rj)d5 


Evaluation of the polynomials in the above integrals leads to the foi* 
lowing forcing boundary conditions for the quadratic element in terms 
of a velocity source which form the right hand side of Eq. (B~3) 


6 : = 


X Cf,(rf(^) + r,rj(5i) + rj(^)) 




«j 


X p-) + r.r.(- To?> * j <’ 5T’ > 


2 

j‘^ 


+ f-.UU- 2ii) + ♦ rj-tTif-)' 


> J 


j M05‘ 


and 


i|(r.-r,)(c^-i?) X [f, • + r.r, + r‘ {- -ry^) ] 


. 2 , 1 


J I 


I ^ 'i '70 


I J 


' 210 ' 




^ V'■^■ 525> * * ^-0* 




B. 5.3 Evaluation of a Pressure Source Fan Plane Condition for a 3~Node 
Triangle (See Figure B~M 


The boundary integral to be evaluated is given by Eq. (B-53) . The 



boundary condition in terms of a radially varying pressure source is 


P • i(r) 
and 

P “ g(r) 

One could solve for i /% in this case) in terms of p from 

n n 2 z 

Eq. (8) only i f 0 to obtain 

- ^ _ [4 -t- p/p] 

and 

; , ~ p/R.) 

^ ♦z 

Substituting Eqs. (13) and (B-67) into Eq. {B~59) one obtains, 

3 ■ - /N r^(c^-;^) ds 

" “ " 

Following the procedure indicated in Section B.5*l and using the 
interpolation functions given by Eq. (B"26) one may write 

9-0 t)9j + C9j 

r = (1 -C)r. + 5 r^ 

^ = 0 'C)^ j 


03) 

( 1 ^) 

(B-67) 

(B-68) 

(B-69) 


(B-70) 


Substituting Eqs. (B- 70 ) into Eq. (B-69) yields 


6 


I 


j 



0-E) 

U) 


/-2 : 2 x 

2 ^ ” *^2 

X [{(l-C)«’, +Cr.}^ : X 

{ta»[(I-E)ij 0jl 

+ r in -E)g| cijl ) 1 (r.-r.)dc 


Evaluating the various polynomials in the above integrals yield the 
following forcing boundary conditions for the linear element in terms 
of a pressure source which form the right hand side of Eq. (6-3) (note 
that the terms involving and 4>j are transported to the left hand 
side of Eq. (B- 3 ) as they are unknown yet) 


S . 


(c^ - 6 


(r.-r.) X 
J ' 


[(r^ + 4j)0'f(-^) + + rj(^)} 


2 , 1 


4j){r^(^) + r.r.(ji) + r?(^)}] 


(B- 71 ) 



B.5.^ Evaluation of a Pressure Source Fan Plane Condition for a 
6~Node Trl^ngi<r (See Figure B-4) 


The boundary Integral to be evaluated for the case of a pressure 
source Is given by Eq. (B~69). Incorporating the Interpolation func~ 
tlwis for quadratic elements given by Eq. (B-30) one may write 

i - (2C^ - 3C + Dij + C)9j 

( 25 ^ - 35 + 0^1 (25^ - C)lj + 

(B-73) 

r ■ (I - c)r. + CTj 
ds “ (r. - r. )d5 

j I 

Substituting Eqs. (B-73) Into Eq. (B-69) yields 


8 . 

I 

j 

I 


_ r ’ ( 25 ^- 35 + 1 ) ^ 

£-0 (2£^-£) 
M£-€^) 


{(l-5)r. <-£r.}' 


/*2 t2. 

(c - 


wl (25^-35+1 )^j + (2C^-£)^j+^«(£-5^)0j] 


- [(25^-3£+l)g.+ (2£^-£)g.+l»(£-r.^)g ] 

P I J ^ 


X (fj - r. )d£ 


Evaluating the various polynomials In the above Integrals yield the 
following forcing boundary conditions for the quadratic element in 
terms of a pressure source which form the right hand side of Eq. (B“3) 
(note that the terms involving i. and are transported to the 

I J 

left hand side of Eq. (B-3) as they are unknown yet) 
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Cc^-;f)(r - r,) 


(r^ + u}|Hr|(|5j) + * '■j*2lo'* 


2. 1 


4 j){r.{- -gi;) + r.r.(- j^) ^ (- ^)) 


*'4j)(rj<75) + * ']<■- 555 )) 

- i^)Uj-r ) 

^ W 


+ 4 ,)fr?(- 55) + r,r.(- -fji) + (- j{) ) 


j+ (j ♦ 4.)(r,(yf3) + 




(c^ - ih (r - r ) 

-k £ i !_ X 


(|i + * V' 


* (r^ + 4jHr,(- -jji) + r,r.(yj|) ♦ r?(yi)) 


2, 1 


l> . .2, 1 


♦<.(—+ u.«, Hr. (-^) + r,rj(.^) + r . (-^55) ) 
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B.6 Boundary Integral Evaluation for Inlet Entrance Plane 
The boundary condition at the Inlet entrance plane In terms of a 
general radially varying radiation Impedance Is expressed In Eqs. (34) 
and (35): 

|^(c6^ - - c x^ = 0 )^ + (341 

c x^ + i^(c (35) 


B.6. 1 Evaluation of Entrance Plane Radiation Condition for a 3~Node 
Triangle (See Figure B~5) 

Noting that along the curve C a is equal to 270“ the boundary 
integral In Eq. (B-6) reduces to 




(B-77) 


Noting that for the configuration In Figure B-5 S and n are in the 
opposite directions to r and z respectively 


B = / N r^[-2$ i ^ + (c^ - ^^)$ ]dr 


(B-78) 


The boundary conditions given by Eqs. (34) and (35) can be solved for 
i and i and on substitution into Eq. (B-78) yields 


r2 l2v 

,2 , 

= ^r N ^ ‘l*z)(‘^‘i> + '^r'^r^ 


Ae 


+ c (-to| + )(]dr (B-79) 


Ae - (c B - 5 + (c X )^ 

e z e 


where 
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Incorporating the interpolation functions given by Eq. (B-26) one may 
wr i te 


<<> * 0 -C)4>. + (B-80) 

r = (1 -C)r. + CTj 

Substituting Eqs. (B-80) into Eq. (B-79) one obtains 


I £=0 ■' 

J 

- $j)/{r. - r^. ) 

- . . 

j. ^ /(j. 0^ _ (|i^ ) (u){ } 


+ 


Tv^T 



♦j)) 


X 



+ (c X {-(*){ (1 - ) 

e I j 




r.)dC 


(B-81) 


Evaluating the various polynomial integrals assuming that the mean flow 
variables and the radiation impedance to be constant over the element 


leads to 
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(c^- 

— 4 e“ X [(ie^-*^)^(r,-r,) (♦ilrf(i) + ^ , r j (-ji) * r] (-j^) 1 


+ r (i I - i j H r? ({) + r . r . (;^) -f r j (fi) ) 


xx^»(r,-r.) {i.[r 2 (^)*r.r.(-fi) + rj(3i)l 

4 jlr^( 2 i) + r.r.(Tj) + r^(^)]) 


+cXg J(S- 82 ) 


Bj = + r.r.(j) -f r^^(|)] 


(- 2 _ 

ra 9 ^-*^)«.(r.- r.) {i.(r2(5i) + r.r.(^)+r](^)] 


+ (^ 0^- - ijUrf (-rj)+ 


-(c x^o.(r,- r.) {i|[r 2 (^) + r.r.(fi) + rj^{^)l 
+c ^ (i .-♦.){ (fi) + r . rj (^) + r j (i) )] (B- 83 ) 


15 ^ 


Since the acoustic potentials are unknown yet, 6j and Bj are transported 
to the left hand sides of the corresponding equations in (B-3) for nodes 
i and j respectively. 

B.6.2 Evaluation of Entrance Plane Radiation Condition for a 6-Node 
Triangle (See Figure B-S) 

The boundary integral to be evaluated is given by Eq. (B-79). 
Incorporating the interpolation functions given by Eqs. (B-30) one may 
wri te 

0 = - 3C + 1)4>. + (2E^ - 

I j *- 

r = (1 - E)r. + ^r. (B-84) 

Substituting Eqs. (B-8I4) into Eq. (B-79) one obtains 


8. . / ' (h\-35*1) , 

\ C=0 (2C - 0 X 

^ ‘iU- 


(c - 
Ae 


{ (c 0 - 5 ) (wi + 

e ^ r r 


+c X (-0)4) + ^ ^ ) } 
e r r 


2 , 

r dr 


- 2 $ 

r Z r 


Evaluating the various polynomials in the above integrals yields 


B. = 

I 


»*^(r - r ) (c^- 

* J ^ 

Ae 


I J 


•J * *j(- 35) * ijC-ij)) 


■t' 
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■ 2 r^ |) + ij(- -^) + 


(B- 


where r « (tj + r^)/Z 


Sj.- 


- 2 . .^2 l 2 x 

r (r. - r. ) (c - (f>^) 
Ae 


(c jw[^.(- -^) + ij (y|) + 


- TF^ ^ ^ 


■c xju)[^.(- 3 ^) + i.(-j|) + 


■" Tr.^r“y 


■ 2 r^ ‘‘‘ "'’ 


(B-86) 


and 




■ 2 r \/-2 l 2 , 

r (r. - r^) (c - <p^) 
Ae 


fc - 1,^) )»(ii(-ii) + ij(f|) + ♦,(-j|)l 


■ (r. -r.) f’ * ♦j'f’ ■" 


I J 


■; |»[+,(-jj) + +j(-i^) + 


■" Tr.-'r.'l '♦i’" I’ * Vf* * i,<“)l|_ 


• J 

2'r^ 4>/4>2 [^j(- |) + ^j(|) + ^JO)] 


(B-87) 
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Since the acoustic potentials are unknown yet, S.» S. and 6 are trans 

I J X 

ported to the left hand sides of the corresponding equations in (B-3) 
for nodes i, J and £ respectively. 


Figure B-3. Linear and Quadratic Elements on Inlet Wall 



Figure B-4. Linear and Quadratic Elements on Fan Plane. 



Figure B-5. Linear and Quadratic Elements on Inlet 
Entranre Plane. 
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APPENDIX C 

ACOUSTIC INTENSITY AND POWER CALCULATIONS 

C. I Derivation of Equation (64) for Acoustic Intensity 
Recall that the assumptions of the theoretical model are: fluid 

is inviscid, non-heat conducting and satisfies the Ideal gas law. The 
fluid flow is isentropic. Hence the steady state stagnation enthalpy 
is constant. The present objective is to evaluate an expression for 
acoustic intensity for a wave propagating in a flow with gradients. 

For any closed surface S fixed in flow the total outward energy 
flow, E is 


E = / N • n ds (C-1) 

S ‘ 

where N is the energy flux vector 

and n is the unit outward normal vector. 

And 


N = J m (C-2) 

where J is the specific stagnation enthalpy 
and m is the mass flux vector. 

Substituting Eq. (C“2) into Eq. (C-1) and taking the time average of 
the equation one obtains 

= < / J rn 

5 




• n ds ’ 


(C-3) 
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I * 

where the operator <( )> ■ Lim •== / { )dt 

T->® -T 

Since S is fitted the time averaging operator and the integral operator 
can be exchanged and noting that n is not a function of time Eq. (C-3) 
reduces to 

<E> = / <Jm > • n dS (C-4) 

S 

The energy flux vector, N, is a second order quantity in terms of 
acoustic perturbations. Hence in the development of an expression for 
N one needs to retain terms up to second order in J and m individually 
rather than considering N as made up of product of only first order per 
turbations. This subtle point has been noted by Zinn^^ in the context 
of evaluating acoustic losses in short rocket motor nozzles. By per- 
forming time averaging operation on the continuity equation containing 
terms up to second order perturbations and also assuming periodic first 
order perturbations, Zinn has shown that the time averaged value of the 
second order perturbation of ip is a constant that can be chosen to be 
equal to zero without any loss of generality. An identical operation 
on Euler's equation for the case of compressible irrotational flows 

shows that the time averaged value of the second order perturbation of 
36 

J is also zero. Hence, the contribution of second order perturba- 
tions of J and m towards the time averaged acoustic intensity is zero. 
Subject to this qualification, one can express J and m as composed of 
only a steady state component and a periodic acoustic component, i .e. , 
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J (x, t) “ J (x) + J' (x,t) 

and (C-5) 

m(x, t) ■ m(x) + m' (x, t) 

Under the isentropic assumption made the steady state stagnation enthalpy 
J is a constant and not a function of the position vector, x. Substl^ 
tuting Eq. (C“5) into Eq. (C-^) and noting that j and m are independent 
of time and that the time averaged values of the periodic acoustic 
quantities are equal to zero, one obtains 

<E> *3 / m*ndS+ / <j'm‘> • n dS (C-6) 

s’ s’ 

The continuity equation corresponding to the steady state quantities 
in the integral form reads 

/ m • n dS = 0 (C-7) 

S ’ 

Hence the first term in Eq. (C-6) contributes nothing towards the out- 
ward energy flow, that is 

<E> * / <J'm'> • n dS (C-8) 

S 

Since J m corresponds to intensity in the absence of acoustic pertur- 
bations one may define the time averaged acoustic intensity as the 
difference in intensities <J m > and J m, that is 

<1 ^. > = <J m > - j m = <J'm*> 

-acoustic 


(C-9) 



}6I 


The specific stagnation enthalpy J is given by 

J - h + Y |vj^ (C-IO) 

where h is the specific enthalpy 
1 2 

and "2 lyi *s specific kinetic energy. 

An acoustic perturbation of Eq. (C-IO) yields 

J' ■ h' + y • y' (C-1 1 ) 

A similar acoustic perturbation of mass flux vector m yields 

m‘ * p y ’ + p ' y (C“ 12) 

The classical Maxwell's relations^^ for the isentropic situation 
yield 

h. » (ill )P’ 

3P s 

P 

and (C-13) 

p. « (l£.)P' . BL 

^ ^9P s -2 

where s is the specific entropy of the system. 

Substituting Eqs. (C-ll) - (C-13) into Eq. (C-9) yields the desired 
expression for the time averaged acoustic intensity 

. 2 - 

P y - - - n' - 

<1 ..>»<p'v'+ r + p (v • v' )v' + (v • v‘ )-^ V > (C- ll») 

-acoustic - --2 ^ -- - - -2 - 

pc 

For the case of an irrotationai and uniform entropy flow, the mean and 
acoustic velocities v and y' in Eq. (C-ll«) can be expressed as gradients 
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of mean and acoustic potentials ^ and 0' respectively. Hence the 
acoustic intensity in the direction of a unit vector a is 


<1 ^ 5 > 

-acoustic 


•a <p'V^'> 


+ p[v^ 


70 • a 2 

a rrj- <p‘ > 

P c 

{V0 • a) 

<70'70'>] ’ a + ■ " ■' " 70 • <p'70'> (64) 


which is the desired relationship. 

The time averaged energy flow across a surface S Is then 


<E> “ / 


<1 . . > 
-acoust I c 


n dS 


(65) 


C.2 Evaluation of Acoustic Power at the inlet Boundaries 

and dB Calculations 

Equation ( 65 ) has to be evaluated at the fan plane. Inlet entrance 
plane and the inlet walls to determine the acoustic power input Into 
the inlet, the acoustic power output of the Inlet and the acoustic p<>rer 
absorbed by the lined inlet walls respectively. 

In this context, it can be shown easily that for two sinusoidal 
functions given by 

q « Re(Q e 
and 

r - Re[R 


where Re[ ] = Real part of [ ], 
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■ Q and R are complex quantities given by 
Q - (i ♦ iii and R - R + ift 
the time averaged value of qr is 

<qr> - Re[Q "] (C-15) 

where R is the complex conjugate of R. 

This result is used repeatedly in evaluating the various time averages 
occurring in Eq. (6^). 

C.2.1 Acoustic Power Input into the Inlet at the Fan Plane 

At the fan plane the unit outward normal points in the positive 
Z-direction. Hence the acoustic power input is 



Equation (C'!?) is evaluated for the finite elements located on the 
fan plane by calculating the time averaged acoustic intensity in terms 
of the element centroidal values of the acoustic variables for the sake 
of simplicity. Hence the acoustic intensity as calculated by Eq. 

(C~l8) over each such element is constant and can be taken out of the 
integral in Eq. (C-l?)* Hence the contribution of one element located 
on the fan plane to Eq. (C-17) is (see Figure 


<E*>. . » 2tt <1® 

input acous 


• ^ / 
tic 2 ■' 


r . 
J 


r dr 


r*r . 
1 


(C-19) 


IT <1® > (r? - r?) 

acoustic z j 1 


Sunwiatlon of Eq. (C"I9) for all the elements located on the fan plane 
yields the total time averaged acoustic power input, i.e., 


<E>. 

input 


■'F 

I 

e=l 


<E®> 


input 


(C-20) 


where N^. Is the number of elements located on the fan plane. 

C.2.2 Acoustic Power Output at the Inlet Entrance Plane 

At the inlet entrance plane the unit outward normal points in 
the negative Z-dIrection. Hence the acoustic power output is 


<E> 


output 


- / 


Entrance plane 


<1 >,dS 

acoustic 2 


(C-21) 


Again Eq. (C-21) can be integrated with respect to e to yield 
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r 

oe 

<E> ^ ^ » -2tt / <1 > rdr 

output ■' acoustic z 

le 


(C-22) 


where r. and r are the Inner and outer radii at the inlet entrance 
I e oe 

plane respectively. 

Equation (C-22) is evaluated for the finite elements located on the 
entrance plane by calculating the time averaged acoustic intensity as 
in Section C.2.1. Hence the contribution of one element located on the 
entrance plane to Eq. (C-22) is (see Figure B-5) 


<E®> 


e . ' 

. . = -2tt <1 > / rdr 

output acoustic z ■' 

j 

= -TT <1® > (r? - r?) 

acoustic z I j 


(C-23) 


Summing of Eq. (C-23) for all elements located on the entrance plane 
yields the total time averaged acoustic pc^er output, i.e.. 


<E> = y <E®> 

output ^_j output 


(C-24) 


where Is the numbei' of elements located on the inlet entrance plane. 
C.2.3 Acoustic Power Absorbed at the Inlet Walls 

The acoustic power absorbed by the inlet walls is 


<E> 

absorbed 


^ r 

J 

Inlet wal 1 


<1 ' dS 

acoustic n 


(C-25) 


S • nee the integrand is independent of 6, Eq. (C-25) can be integrated 
with respect to 6 to yield (see Figure B*3) 
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<E> . u j 2 tt / 
absorbed ■' 

s 


<1 > rds 

acoustic n 


(C-26) 


where s Is the natural coordinate along the inlet wall. Substituting 


Eq. Into Eq. (C-26) leads to the expression for<E>_ 


absorbed 


for a single element located on the inlet wal 


<E®> 


absorbed 


= 2TI J 


4=0 


<1 .. > [ (l-4)r.+4r.]L . .d^ 

acoustic n i ^ j i j ^ 

(C-27) 


Noting that at the Inlet walls the normal component of the mean flow 
velocity is zero, Eq. (64) reduces to the following expression for the 
acoustic intensity normal to the inlet wall 


<1 . . > 

acoustic n 


= <p' 


3n 


p(v^ ■ V(!)' ) 


iil 

9n 


(C-28) 


Evaluating the vector dot product in terms of the natural coordinates 
(s'n) and using Eq. (C-l5)f Eq. (C-28) reduces to 


I =:^[(p^+p<i)) + 

acoustic n 2 ^n n 


^s^n^^ (C-29) 


Equation (C-27) is evaluated for the finite elements located on the 
inlet wall by calculating the normal component of the time averaged 
acoustic intensity in terms of the centroidal values of the acoustic 
variables using Eq. (C-29). Hence the contribution of one element 
located on the inlet wall to Eq. (C-27) is 


rS , e 

*^E > . u j = 2 tt< I . . > 
absorbed acoustic n 



= TT 


acoustic n 


L . . ( r . + r . ) 

ij ' J 


(C-30) 
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Summing of Eq. (C-30) for all elements located on the inlet wall yields 
the total time averaged acoustic power absorbed by the Inlet wall, i.e., 


N 

0) ^ 

^^^absorbed " ^ ^absorbed 

e=l 


(C-31) 


where N is the number of elements located on the upper and lower walls 
of the Inlet. An obvious conclusion from Eq. (C-29) is, the normal 
component of the acoustic intensity at a hard wall is zero since * 

n 

and are prescribed to 'e zero there. Hence for a hard wall the 
acoustic power absorbed is zero. 

C.2.4 DB Calculations and the Principle of Conservation of Acoustic 

Enerai 

Once the energy fluxes at the inlet boundaries are determined by 

Eqs. (C-20), (C-24) and (C“31) the dB calculations are performed to 

estimate the effectiveness of a liner. The dB , for a liner is 

reduction 

defined by 

dB . _ = 10 log.. 

reduction lO 

The orinciple of conservation of acoustic energy states that the 
acoustic power input should equal the sum of acoustic power output and 
the acoustic power absorbed, i.e.. 


<E> . 


<E> 


nput 


output 


( 66 ) 


I nput 


<E> ^ + 

output 


<E> 


absorbed 


(C-32) 


For a hard walled inlet Eq. (C~32) reduces to 
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<E>, 


input 


<E> 

output 


and hence 


dB 


reduction 


Hard walled 
Inlet 


= 0 


(C-33) 
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APPENDIX 0 

GEOMETRICAL AND MEAN FLOW DATA FOR THE QCSEE AND 
BELLMOUTH INLETS 

The quadratic triangul ization scheme Is used for predicting the 
attenuation of sound due to liners on the inlets' upper walls (see 
Figures 4a and 4b). As one may observe the total number of nodes 
equals 289 and the total number of elements equals 126 for both the 
QCSEE and Bellmouth inlet triangul ization schemes. The total number 
of corner nodes equals 82 and the total number of mid-side nodes equals 
207. The geometry of the inlets is specified in terms of the corner 
node number, IL and its Z- and r- coordinates, Z()L) and R(IL) respec- 
tively. The "compressible" two dimensional axisymmetric mean flow at 
the corner node number, IL in terms of the axial and radial velocity 
components, (}i^(lL) and (f>^(lL) correspond to an average exit Mach number 
of 0.52. The free stream Mach number, M„ equals 0.12 for the QCSEE 


inlet and M is 0.0 for the Bellmouth inlet. 
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Table 0-1. 

Geometrical and 
the QCSEE Inlet 

Mean Flcn^t Data 

for 


(M « 0. 12, 

w 

— 0.52{ See 

Figure 4-a) 


'l 

Z(1L) 

R(1L) 



1 

0.0 

1 . 0500 

0.0000 

-0.5387 

3 

0.0 

0.7875 

0.4364 

-0.2204 

5 

0.0 

0.5250 

0.4401 

-0.0958 

7 

0.0 

0.2625 

0.4330 

-0.0276 

9 

0.0 

0.0000 

0.4315 

0.0000 

20 

0.1 

0.9182 

0.6271 

-0.3934 

22 

0.1 

0.8034 

0.5993 

-0.3026 

2 k 

0.1 

0.5739 

0.5166 

-0.1129 

26 

0.1 

0.3443 

0.4850 

-0.0415 

28 

0. 1 

0.1148 

0.4766 

0.0057 

30 

0. 1 

0. 0000 

0.4752 

0.0000 

k ] 

0.2 

0.8768 

0.7618 

-0.2545 

k 3 

0.2 

0.6576 

0.6283 

-0.1264 

kS 

0.2 

0.4384 

0.5484 

-0.0577 

kl 

0.2 

0. 2192 

0.5232 

-0.001 1 

kS 

0.2 

0.0000 

0.5146 

0.0000 

60 

0.3 

0.8564 

0.8215 

-0.1358 

62 

0.3 

0.7494 

0.7196 

-0.0868 

(>k 

0.3 

0.5353 

0.6047 

-0.0502 

66 

0.3 

0.3212 

0.5646 

-0.0101 

68 

0.3 

0.1071 

0.5600 

0.0400 

70 

0.3 

0.0000 

0.5590 

0.0000 

81 

0.4 

0. 8500 

0.7894 

-0.0205 

83 

0.4 

0.6375 

0.6930 

-0.0353 

85 

0.4 

0.4250 

0.6146 

-0.0234 

87 

0.4 

0. 21?5 

0.5845 

0.0115 

89 

0.4 

0. 0000 

0.5759 

0.0000 
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Table O-l. (Continued) 


'l 

Z(IL) 

R(IL) 



100 

0.5 

0.8516 

0.7395 

0.0188 

102 

0.5 

0.7452 

0.7059 

0.0037 

] 0 k 

0.5 

0.5323 

0.6355 

-0.0086 

106 

0.5 

0.3194 

0.6131 

-0.0042 

lod 

0.5 

0. 1065 

0.5948 

0.0505 

110 

0.5 

0. 0000 

0.5942 

0.0000 

121 

0.7 

0.8638 

0.6648 

0.0550 

123 

0.7 

0.6479 

0.6340 

0.0291 

125 

0.7 

0.6319 

0.6067 

0.0185 

127 

0.7 

0.2160 

0.5892 

0.0196 

129 

0.7 

0.0000 

0 5865 

0.0000 

140 

0.9 

0.8847 

0.5969 

0.0710 

11«2 

0.9 

0.7741 

0.5970 

O.O6O3 


0.9 

0.5529 

0.5746 

0.0408 

Ili6 

0.9 

0.3318 

0.5621 

0.0346 

148 

0.9 

0. 1 106 

0.5472 

0.0625 

150 

0.9 

0.0000 

0.5449 

0.0000 

161 

1.1 

0.9109 

0.5401 

0.0751 

163 

1.1 

0.6832 

0.5425 

0.0616 

165 

1.1 

0.6555 

0.5246 

0.0559 

167 

1.1 

0.2277 

0.4953 

0.0491 

169 

1.1 

0.0000 

0.4608 

0.0000 

180 

1.3 

0.9389 

0.5023 

0.0699 

182 

1.3 

0.8215 

0.5131 

0.0680 

184 

1.3 

0.5868 

0.5005 

0.0773 

186 

1.3 

0.3521 

0.4734 

0.0900 

188 

1.3 

0. 1174 

0.3461 

0. 1180 

190 

1.3 

0.0000 

0.2284 

0.0000 
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Table D-:I. (Continued) 


'l 

2(IL) 

R(IL) 



201 

1.4 

0.9525 

0.4876 

0.0647 

203 

1.4 

0. 7144 

0.5042 

0.0719 

205 

1.4 

0.4763 

0.4880 

0.0988 

207 

1.4 

0.2381 

0.3829 

0.1657 

209 

1.4 

0.0000 

0.0000 

0.0000 

220 

lo 

0. 9652 

0.4817 

0.0593 

222 

1.5 

0.8722 

0.4985 

0.0595 

224 

1.5 

0.6862 

0.5067 

0.0691 

226 

1.5 

0.5002 

0.5036 

0.0981 

228 

1.5 

0.3142 

0.4433 

0.2121 

230 

1.5 

0.2212 

0.3765 

0.3361 

241 

1.6 

0.9765 

0.4748 

0.0524 

243 

1.6 

0.8069 

0.4999 

0.0540 

245 

1.6 

0.6373 

0.5142 

0.0763 

247 

1.6 

0.4676 

0.5192 

0. 1474 

249 

1.6 

0.2980 

0.4970 

0.2980 

260 

1.8 

0.9935 

0.4750 

0.03)6 

262 

1.8 

0.9165 

0.4853 

0.0305 

264 

1.8 

0.7624 

0.5127 

0.0391 

266 

1.8 

O.6O83 

0.5408 

0.0626 

268 

1.8 

0.4542 

0.5977 

0. 1282 

270 

1.8 

0.3771 

0.6177 

0.1583 

281 

2.0 

1.000 

0.5053 

0.0000 

283 

2.0 

0.8500 

0.5124 

0.0108 

285 

2.0 

0.7000 

0.5316 

0.0205 

287 

2.0 

0.5500 

0.5634 

0.0251 

289 

2.0 

0.4000 

0.5887 

0.0000 
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Table D-2. Geometrical and Mean Flow Data for 
the Bel Imouth Inlet 

(M » 0.0, M “ 0.52; See Figure ^-b) 

00 Q 


'l 

Z(IL) 

R(IL) 



1 

0.0000 

1 . 0000 

0.0000 

-0.2020 

3 

0.0000 

0.7500 

0.2239 

-0.1860 

5 

0. 0000 

0.5000 

0.2460 

-0.0885 

7 

0. 0000 

0.2500 

0.2459 

-0.0381 

9 

0. 0000 

0.0000 

0.2457 

0.0000 

20 

0.0690 

0.8147 

0.1770 

-0.2932 

22 

0.0690 

0.7129 

0.2891 

-0.1684 

Ik 

0.0690 

0.5092 

0.2764 

-0.0898 

26 

0.0690 

0.3055 

0.2690 

-0.0466 

28 

0.0690 

0. 1018 

0.2655 

-0.0146 

30 

0.0690 

0.0000 

0.2620 

0.0000 

41 

0. 1379 

0.7716 

0.3492 

-0. 1856 

43 

0.1379 

0.5787 

0.3144 

-0.1033 

45 

0.1379 

0.3858 

0.2956 

-0.0581 

47 

0.1379 

0. 1929 

0.2868 

-0.0265 

49 

0. 1379 

0.0000 

0.2781 

0.0000 

60 

0.2759 

0.7198 

0.3908 

-0.0978 

62 

0.2759 

0.6298 

0.3820 

-0.0856 

64 

0.2759 

0.4499 

0.3452 

-0.0546 

66 

0.2759 

0.2699 

0.3267 

-0.0306 

68 

0.2759 

0.0900 

0.3191 

-0.0099 

70 

0.2759 

0.0000 

0.3115 

0.0000 

81 

0.4138 

0.6897 

0.4306 

-0.0537 

83 

0.4138 

0.5173 

0.3853 

-0.0361 

85 

0.4138 

0.3449 

0.3597 

-0.0261 

87 

0.4138 

0. 1724 

0.3478 

-0.0134 

89 

0.4138 

0.0000 

0.3358 

0.0000 


17 ^ 


Table D*2. (Continued) 


'l 

Z(IL) 

R(IL) 


i^(lL) 

too 

0.5517 

0.6897 

0.3818 

0.0000 

102 

0.5517 

0.6035 

0.3984 

-0.0057 

104 

0.5517 

0.431 1 

0.3808 

-0.0150 

106 

0.5517 

0.2586 

0.3682 

-0.0113 

108 

0.5517 

0.0862 

0.3621 

-0.0040 

110 

0.5517 

0.0000 

0.3559 

0.0000 

121 

0.8276 

0.6897 

0.3640 

0.0000 

123 

0.8276 

0.5173 

0.3858 

-0.001 1 

125 

0.8276 

0.3449 

0.3811 

-0.0011 

127 

0.8276 

0. 1724 

0.3762 

-0.0005 

129 

0.8276 

0.0000 

0.3714 

0.0000 

140 

1.1034 

0.6897 

0.3661 

0.0000 

142 

1.1034 

0.6035 

0,3903 

0.0018 

144 

1.1034 

0.4311 

0.3852 

0.0075 

146 

1.1034 

0.2586 

0.3743 

0.0110 

148 

1.1034 

0,0862 

0.3618 

0.0061 

150 

1.1034 

0.0000 

0.3492 

0.0000 

161 

1.2414 

0.6897 

0.3810 

0.0000 

163 

1.2414 

0.5173 

0.3942 

0.0090 

165 

1.2414 

0.3449 

0.3803 

0.0225 

167 

1.2414 

0. 1724 

0.3465 

0.0295 

169 

1.2414 

0.0000 

0.3167 

0.0000 

180 

1.3793 

0.6897 

0.3986 

0.0000 

182 

1.3793 

0.6035 

0.4089 

0.0075 

184 

1.3793 

0.4311 

0.3976 

0.0290 

186 

1 . 3793 

0.2586 

0.3633 

0.0655 

188 

'.3793 

0.0862 

0.2123 

0. 1121 

190 

1.3793 

0. 0000 

0.0000 

0.0000 
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Table D-2. (Continued) 


'l 

Z(IL) 

R(1L) 


0^(10 

201 

1 . 4483 

0.6897 

0.4046 

0.0000 

203 

1 . 4483 

0.5485 

0.4134 

0.0168 

205 

1 . 4483 

0.4074 

0.4020 

0.0413 

207 

1 . 4483 

0.2662 

0.3731 

0.0867 

209 

1.4483 

0.1250 

0. 1880 

0.2350 

220 

1.5172 

0.6897 

0.4121 

0.0000 

222 

1.5172 

0.6250 

0.4240 

0.0084 

22 k 

1.5172 

0.4957 

0.4188 

0.0294 

226 

1.5172 

0.3664 

0.4060 

0.0616 

228 

1.5172 

0.2371 

0.3794 

0.1282 

230 

1.5172 

0.1724 

0.3480 

0.2067 

2i«1 

1.6552 

0.6897 

0.4304 

0.0000 

243 

1.6552 

0.5755 

0.4426 

0.0194 

245 

1.6552 

0.4613 

0.4385 

0.0460 

247 

1.6552 

0.3470 

0.4283 

0.0871 

249 

1.6552 

0.2328 

0.4107 

0.1540 

260 

1.7931 

0.6897 

0.4514 

0.0000 

262 

1.7931 

0.6897 

0.4514 

0.0000 

264 

1.7931 

0.5378 

0.465’7 

O.O2C9 

266 

1.7931 

0.4365 

0.4664 

0.0589 

268 

1.7931 

0.3352 

0.4656 

0. 1068 

270 

1.7931 

0.2845 

0.4587 

0.1434 

281 

2.0000 

0.6897 

0.4745 

0,0000 

283 

2,0000 

0.6013 

0.4940 

0.0110 

285 

2.0000 

0.5130 

0.5006 

0.0213 

287 

2.0000 

0.4246 

0.5183 

0.0352 

289 

2.0000 

0.3362 

0.55’3 

0.0421 
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